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ABSTRACT 


A simplified  two~degr ee-of-freedom  simulation  in 
heave  and  pitch  of  the  XR-3  Captured  Air  Bubble  (CAB) 
testcraft  is  developed  first  in  a non-linear  form  and 
is  subsequently  linearized  about  the  steady-state 
operating  point.  The  model  is  validated  against  the 
six-degree-of-f reedcm  XR-3  Loads  and  Motions  program 
in  the  time  domain;  the  linear  system  is  then 
transformed  into  the  frequency  domain  by  complex 
matrix  inversion.  A comparison  is  made  between  the 
linear  system  frequency  response  predicted  by  the 
two-degree-of- freedom  model  and  that  of  the  Loads  and 
Motions  program,  a-s  well  as  actual  testcraft  data, 
reduced  to  the  frequency  domain  by  a Fast  Fourier 
Transform  (FFT)  technique.  This  comparison  of 
frequency  response  curves  highlights  a non-linear  mode 
of  CAB  craft  behavior  of  possible  interest  to 
habitability  and  seaworthiness  studies. 
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(lof) 


Dpsmb 


Partial  derivative  of  Bowseal  pitch 
moment  with  respect  to  db 

( ft-lbf/lbm) 


Dpsth 1 


Partial  derivative  of  Bowseal  pitch 
moment  with  resoect  to  Tneta 

( ft-lof) 


Dpsth2 


Partial  derivative  of  Sternseal 
pitch  moment  with  respect  to  Theta 

(f  t*lbf> 


Dpsz  1 


Partial  derivative  oi  Bowseal  pitch 
moment  with  respect  to  z 

( Ibf) 
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30.  Dps22 


Partial  derivative  of  Sternseal 
pitch  moment  with  respect  to  z 

(lof) 


31.  D£ab 


Partial  derivative  of  pitch 
acceleration  with  respect  to  Mb  ^ 

( 1/(lbm-sec)) 


32.  D£th 


33 . Dpz 


34.  Dzmb 


35.  Dzth 


36 . 022 


Partial  derivative  of 
acceleration  with  respect  to 


pitch 

Ineta 

-2 

( sec) 


Partial  derivative 
acceleration  with  respec 

( 


of  pitch 
t to  z 

2 

1/  (ft-sec)  ) 


Partial  derivative  of  vertical 
acceleration  with  respect  to  Mb 

2 

( ft/lbifl-sec) 


Partial  derivative  of  of  vertical 
acceleration  with  respect  to  Iheta^ 

(ft/sec) 


Partial  derivative  of  vertical 
acceleration  with  respect  to  z 

-2 

( sec) 


I 
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INTRODuCTIUN 


Ih>?  rr34U«?ncy  aomaia  aroviJes  a asaiul  Ti-''-tao'.l  of  atiiy 
for  sy:jt£:n;i  such  2o  tlie  Cdtitarod  aii:  Buucla  (CAa)  tyae 
surtdc.2  eiiiiCT.  ship,  wtiicii  ara  axciced  by  torsos  «noj*'- 
statistics  of  asuiitude  a:rl  tre^aoiicy  j.ay  o'-’  wail- tnow  r: , ;t 
wnosc  title  doHaia  piidsi  roi  It  ionsni -5  ir-j  ta.co 
indetertiir.aut.  In  audition,  tno  troi'joncy  r-sv-ous?  or  suca 
niga-speed  ships  sheuid  ae  closely  exaniiicd  for  rosor,  mess 
near  wuve  encounter  frequencies  most  iiKoiy  tc  jo 
exq;  or  ien  ced . Due  to  its  reLativoly  shalicw  uri*.'t,  tho  C,\J 
crart  ca;i  oasia/  be  -arcitoi  into  <i  no;.-iin<;ar  r-rjicn  ■''f 
operation  in  which  slii'Wiil  ^apjin.j  or  -a ve/o io i. u a up;j-\L 
surtacj  contact  may  occur;  ti:oss  non-iineat  itios  j i y ‘xcit^'* 
modes!  of  operation  which  are  of  irioortance  to  na oi  ta  J ii i * y , 
and  may  limit  the  craft's  s;;eod  during  sea-stato  operation. 


Eootn  in  ref£lj  analyzed  the  frjauei'.cy  uomain  response 
of  the  Xa-3  LAi5  craft  utilizinij  the  i.oaus  and  .'Cotton  pLo;ri::i 
develoqjed  by  uce  an  ics,  I nc . , for  the  j^-li  1JJJ  c;u  in' 
medified  extensively  jy  Leo  and  doncai  (ref 

(rerfjJ)  , ror  the  Xh-3.  how-ver,  ti:o  Loans  a n.i  cotian. 
program  is  an  '^xtremeiy  inefficient  tool  to  itiiizr  f ir 
r’ourier  series  analysis;  not  only  ar  > «avss  u u re  r ica  1 1 y 
integrated  alcng  the  hull  utilizing  a taol^  ioox-^.,  one.' 
for  eacn  frequency  coiiq'onent,  ani  four  tires  .lor  eacr.  ti.i.-_- 
step,  Jut  the  higner  freguency  comnonen  ts  also  iriv-  t.ie 
variable-steq  integrators  to  extre.3''ly  sa-.ill  stt-p  size-,. 
The  nat  eftect  ot  tnis  tecnnipi*  r.i  tn-  -.corfitant 
utilization  or  (.P'J  time  in  a co-ru'-  t si:, uiation  or 
sea-stat  ’ conditions.  in  a i lit  ion,  r..-  ‘ir,^  constant;  oi 
the  crane  nesponsc-  are  r 'iativ--iy  1 ^ n n a to  i*..'.; , 
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particularly  pitch;  therefore  still  further  CP'J  tiae  is 
wasted  attaining  steady-state  conditions.  This  inefficiency 
of  CPU  tiaie  has  been  noted  by  Mitchell  in  ref  fuj  , who 
documented  over  40  percent  of  CPU  time  utilized  oy  the  XK-3 
Loads  and  Motion  program  to  the  WAVES  suoroutine. 

Accordingly,  one  is  reduced  to  performing  a point-by-point 
frequency  analysis  with  a single  wave  frequency  at  a time. 

Since,  under  these  conditions,  the  system  is  approximately 

linear,  and  superposition  is  valid  for  sufficiently  small 

wave  amplitudes,  the  utilization  of  a linearized,  albeit  \ 

highly  simplified,  program  presents  itself  as  a most 

desirable  alternative.  The  effect  of  non-linearities  may  be 

added  to  the  linear  model  with  a describing  function,  and 

the  system  investigated  in  the  complex  plane  utilizing  the 

modified  Nyquist  criteria.  Finally,  a simplified  model 

provides  a starting  point  for  the  investigation  of 

sensitivity  of  craft  operation  to  the  variation  of  assorted 

parameters  through  the  technique  of  root  locus  plotting,  and 

can  be  the  heart  of  the  development  of  any  automatic  control 

system,  optimal  or  classical. 

Gerba  and  Thaler,  in  ref  , initiated  the  development  j 

of  the  linearized  XR-3  model  in  heave  only.  In  this  work, 
that  model  is  further  developed  to  include  the  pitch  mode,  \ 

first  in  a non-linear  form,  from  which  a linearized  version  j 

is  derived.  The  wave  forcing  functions  acting  along  the  ‘j 

hull  of  the  XR-3  are  developed  in  the  form  of  a guasi-linear  j 

B vector  in  state-space,  and  the  linear  system  is  then  | 

transformed  into  the  frequency  domain  by  complex  matrix 
inversion,  with  significant  savings  in  CPU  time.  This 
linear  frequency  response  is  then  compared  with  ref  [<J  >-o 
provide  additional  insights  into  the  XR-3  system  behavior 
not  readily  apparent  in  the  more  complex  Loads  and  Motions  j 

program.  | 

i 

A technique  utilizing  a mini-computer  and  spectral  j 
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analyzer  is  developed  for  the  fast  and  efficient  reduction 
of  recorded  XS-3  data  from  the  time  domain  to  the  frequency 
domain  utilizing  a Fast  Fourier  Transform  (FFT)  . Actual 
craft  data  is  compared  with  the  predicted  linear  system 
frequency  response. 
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II.  NON-LINEAR  XR^S  HSAVS-PITCH  EQUATIONS  OF  MOTION 


This  chapter  develops  the  non-linear  eyaatioas  of 
motions  of  the  XR-3  in  heave  and  pitch,  following  tne 
approach  utilized  by.  Gerba  and  Thaler  in  the  development  of 
the  simplified  heave-only  model  in  ref  CO- 

A.  ASSUdEIICNS 

As  in  ref  fsj , the  following  assumptions  are  made: 

1.  constant  leakage  area; 

2.  equilibrium  speed  conditions; 

3.  plenum  pressure  and  volume  related  through  an  adiabatic 
process ; 

The  following  assumptions  and  limiting  conditions  were 
removed; 


4.  vertically  coincident  geometric  center , center  of 

pressure  and  center  of  gravity; 

5.  negligible  pitch  changes; 

6.  negligible  wave  action; 

7.  negligible  added  mass  effects; 

and  the  foilcwing  assumptions  are  added; 
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negligible  wave  action  within  plenum; 


9.  wave  action  due  to  head-on  or  following  saas 
(negligible  roll  or  yaw); 

In  addition  some  simplifying  assumptions  regarding  the  craft 
geometry  were  made.  Fig  1 shows  the  simplified  geometry 
of  the  craft  assumed  in  this  model.  Note  that  Center  of 
Gravity  (C.G. ) and  Center  of  Pressure  (C.P.)  are  aroitrarily 
located  and  not  to  scale. 

B.  COORDINATE  SXSTEfl 

Fig  1 shows  the  coordinate  system  used  in  this  modal. 
Specifically,  the  coordinate  system  is  attacned  to  the 
geometric  center  of  the  craft  at  the  calm  waterline,  and 
oriented  in  a right-hand  system  such  that  x is  positive 
forward,  y is  positive  to  the  right,  and  z is  positive  down. 
Angles  are  measured  such  that  a roll  to  the  right,  a pitch 
up  and  a yaw  right  are  positive.  Note,  however,  that 
translations  in  the  y direction  and  roll  and  yaw  are  not 
included  in  this  model. 

Draft,  the  distance  of  the  keel  below  the  calm  water 
line,  is  the  distance  of  the  C.G.  above  the  cairn  water  line, 
2,  plus  the  distance  of  the  keel  below  the  C.G.,  Zs,  or: 

Ld  =z+Zs 


(II-1) 


C.  SIDEWALL  GSOilETBY 
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The  constant  sidewall  area  as  used  by  Gerba  and  Thaler 
yields  a good  approximation  to  craft  motion  in  the  heave  and 
pressure  modes.  However  ,it  is  a source  of  strong  coupling 
between  draft  and  pitch,  and  a more  detailed  modeling  is 
necessary  to  achieve  good  agreement  with  the  Loads  and 
Motion  Program,  as  well  as  known  craft  behavior  in  pitch. 

Fig  1 shows  the  cross-sectional  detail  of  assumed 
sidewall  geometry.  The  fore  and  aft  sections  are  identical 
with  the  exception  of  deadrise  angle,  that  angle  being  60° 
at  the  bow  and  80°  at  the  stern.  The  actual  sidewall  width 
used  for  computation  of  buoyancy  forces  is  the  effective 
sidewall  width;  this  is  the  average  value  of  the  waterline 
sidewall  width  and  the  keel-line  width, WslO  or  VJs20  (fore  or 
aft  section) . This  is  one-half  the  draft  averaged  from  the 
C.G.  to  the  bow  and  stern  respectively,  divided  by  the 
tangent  of  the  appropriate  deadrise  angle  plus  the  keel-line 
width  WslO  or  Ws20.  This  value  is  constrained  between  zero 
and  Ws,  the  maximum  sidewall  width.  Thus: 

Ldbarl  = Ld- (L/2-Xcg) *Theta 


(II-2) 


Wsl  =Ldbar  1/ (2*tan  (Dr  1)  ) +Ws  10 

IF  (Wsl  .GT. Ws) Ws1  = Ws 
IF (Hsi .LT. 0. 0) Ws1=0. 0 


and  for  the  stern  section: 


(II-3) 
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Ldbar2 


= Ld+ (L/2+Xcg)  *Theta 


(II-4) 


Ws2  =Ldbar2/(2*tan  (dr2)  ) +«fs20 

IF  (Hs2.cr. Ws) Ws2=Ws 
IF  <Ms2.LT. 0)  Ws2=0 .0 


(II-5) 


D.  SYSTEM  DYNAMICS 


The  equations  of  the  system  dynamics  are  developed  in 
this  section  and  are  an  extension  of  the  model  of  rer  fsj , 
modified  as  necessary;  the  remaining  equations  are  derived 
from  first  principles. 


1 • Pressure  and  M^s  Flo«  Rate 

Plenum  pressure  and  plenum  air  mass  are  calculated 

using  the  adiabatic  Ian.  Using  Ld  as  the  draft  at  the  bow, 

1 

and  Ld  as  the  draft  at  the  stern,  the  average  draft  overall 
2 

becomes: 

Ldbar  =(Ld+Ld)/2 

1 2 

= (Ld-  (L/2-Xcg)  *tan  (Theta) 
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♦ Ld*  (L/2  + Xcg)  *tan  (Theta)  ) /2 


= Ld  + Xcg*tan  (Theta) 


(II-6) 

which,  when  multiplied  by  the  Plenum  area, yields  the 
submerged  Plenum  volume.  Subtracting  this  from  Plenum 
nominal  volume  yields  the  net  Plenum  volume,  or: 


Vb 


= Vn-Ab*  (Ld+Xcg*tan  (theta) ) 


(II-7) 

The  adialatic  relationship  between  pressure  and  volume 
yields: 


Pb 


= Pa*  (Mb/(Vb*Rhoa)  ) 


Gamma 


(II-8) 


and 

Pbbar 


=Pb-Pa 


(II-9) 

Volume  flow  into  the  Plenum  is,  assuming  a constant  fan  map: 
Cin  =N*  (QlO-Pbbar) 


(II«10) 


and  flow  cut  is: 


Qout 


=Cn*Al*  (2*Pbbar/8hoa) 


1/2 


(11-11) 


thus  the  net  mass  flow  rate  within  the  plenum  is: 
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d/  (Mb)  =fihoa>i‘ (Qin-Qout) 
/dt 


(11-12) 


2.  SuD  of  Forces 

Fig  2 shows  the  forces  and  their  respective  monrent 
arms  acting  on  the  XH-3.  The  sum  of  vertical  forces  yields 
the  vertical  C.G.  acceleration: 

2 

d/  2 (z)  i (Hpres+Hbuoy+Hseal+Hdamp+Hplan+W) /M 

/dt 


(11-13) 

The  individual  force  contr ibutions  may  be  calculated 
separately  and  their  total  contributions  summed- 


a.  Pressure  Forces 


The  plenum  pressure  acting  on  the  roof 
plenum  produces  a pressure  force  acting  at  the  C.?.; 


of 


the 


Hpres  =-Ab*pbbar 


(11-14) 


t.  Buoyant  Forces 


For  purposes  of  calculating  buoyant  moments  in 
the  subsequent  calculations,  it  is  convenient  to  dissociate 
the  buoyant  force  into  two  components, Hbuoy  1 and  Hbuoyl,  the 
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buoyant  forces  generated  by  the  submerged  volume  fore  and 
aft  of  the  C.G.  Buoyant  force  is  proportional  to  submerged 
volume  multiplied  by  the  density  of  the  fluid  and  the 
gravitational  constant.  This  term  must  be  doubled  to  account 
for  both  sidewalls..  The  average  draft  forward  of  the  C.G. 
is: 

Ldbarl  =(Ld+Ld1)/2 

= Ld- (L/2-Xcg)  *tan (Theta) /2 


(11-15) 


Thus  the  forward  buoyant  force  is; 

Htuoyl  =-2*Rho*G*length*draft*sidewall  width 


= -2  *Rho*G* (L/2-Xcg) *Ldbar1*Ms1 


(11-16) 

The  aft  average  draft  is: 

Ldbar2  =Ld'»- (L/2+Xcg)  *tan  (Theta) /2 

(11-17) 

and  the  aft  buoyant  force  is; 

Hbuoy2  =-2  *Rho*G* (L/2+Xcg) ♦Ldbar2*Ws2 

(11-19) 

and : 

Hbuoy  =Hbuoy 1 ♦Hbuoy2 


(11-19) 


c.  Bow  and  Sternseal  Forces 
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Fig  3 shows  the  detail  of  the  assumed  bowseai  and  sternseal 
geometry  and  water  surface  interactions.  The  seai  is 
assumed  to  trail  from  the  upper  hinge  point  at  a constant 
angle  of  31. 6°  until  it  contacts  the  waterline.  From  that 
point  on  it  is  assumed  to  lie  tangent  to  the  waterline.  The 
length  of  bowseai  in  contact  with  the  water  may  be  seen  to 
be : 

Xseall  =Ld  1/cos  (3  1 . 6°) 

= 1. 65*  (Ld-  (L/2-Xcg)  *tan  (Theta) ) 


(II-20) 

Pbtar  exerts  a force  upon  this  wetted  area  which  is  resisted 
by  an  upward  reaction  from  the  water  surface.  This  reaction 
is  transmitted  to  the  craft  as; 

Hseall  =-Pbbar *Width*Xseal1 


(11-21) 

As  may  he  seen  from  Fig  3,  the  sternseal 

geometry  is  basically  identical  to  that  of  the  bowseai; 
however,  the  mechanism  for  force  generation  is  different  in 
that  the  sternseal  rides  atop  an  exhausting  jet  of  plenum 
air.  Therefore,  the  pressure  differential  is  not  that  of 
Plenum  pressure  over  atmospheric,  but  over  the  reduced 
pressure  of  the  high-velocity  jet.  This  is  a complicated 

function  of  leakage  area  and  other  factors,  but  may  be 
assumed  approximately  constant.  This  model  assumes 

2 

approximately  2.0  Ibf/ft  differential  pressure  between  the 
sternseal  and  plenum  pressure. 

The  length  of  sternseal  in  contact  witn  the 

water  may  be  seen  to  be: 
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Xseal2 


= Ld2/cos(31.6'>) 


= 1.65*  (Ld+  (L/2  + Xcg)  *taa  (Theta) ) 


(11-22) 

the  sternseal  force  is  then  the  differential  pressure, Sealp, 
tines  the  sternseal  area  ,or: 

Hseal2  =-Sealp*Width*Xseal2 


(11-23) 


d.  Damping  Forces 


Damping  forces  are  developed  for  information 
only,  as  Taylor's  series  expansion  of  these  forces  about  the 
zero-rate  equilibrium  point  causes  these  equations  to  vanish 
in  the  linear  model. 

A body  at  motion  within  a fluid  experiences  a 
frictional  force; 

f =-C*  (fiho/2) *V*1 VJ*S 

f f 

where  C is  a non-dimensional  friction  coefficient  based  upon 

shape,  roughness,  and  Reynold's  Number.  Consistent  with  the 
six-degree-of-f reedom  Loads  and  Motion  program,  C^is  taker. 

to  be  .9.  The  vertical  velocity  is  taken  to  be  the  vertical 

velocity  cf  the  C. G. , 2dot,  plus  the  tangential  velocity  of 
the  angular  velocity,  which  varies  along  the  hull  as; 
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V 00 


=Zdot+  (Xcg-x)  *0taega 


(11-24) 

The  differential  surface  wetted  area  also  is  a function  of 
X,  as  the  depth  of  iamersion  varies  along  the  hull,  so  that: 

dS  =Hs*  (Ld-(x-Xcg)  *tan  (Theta)  ) 

(11-25) 

When  fflulriplied  oy  the  absolute  value  of  velocity,  tiaes 
itself,  times  the  differential  wetted  area,  the  resulting 
expression  may  be  integrated  numerically  along  the  null. 
However,  this  is  most  wasteful  of  CPU  time,  and  the  absolute 
value  term  renders  this  expression  intractable  to  closed 
analytical  solution..  In  addition,  the  term  is  guite  small, 
and  in  fact  must  be,  if  pitch  angles  are  to  remain  within 
the  small  angle  assumptions,  and  no  sidewall  gapping  occur. 
Inerefore,  this  term  is  approximated  by  computing  the 
average  velocity  fore  and  aft  of  the  center  of  gravity,  and 
the  average  wetted  surface  fore  and  aft,  previously  computed 
as  Wsl^Ldbar"!  and  Ws2*Ldbar2.  For  purposes  of  pitch 
calculations  to  be  discussed  subsequently,  these  two  forces 
are  assumed  to  act  at  two-thirds  the  length  fore  and  aft  the 
C.G.;  this  yields  closer  agreement  to  the  numerical 
solution  by  weighting  the  higher  velocities  which  cccur  at 
the  bow  and  stern.  Vbarl  and  Vbar2  may  now  be  defined; 

Wbar  1 =Zdot- (1/2- Xcg) ’•‘Omega /2 

(11-26) 

Vbar2  =Zdot+ (L/2  + Xcg) ’•‘Omega/2 

(11-27) 

Since  there  are  four  surfaces  forward  and  aft  of  the  center 
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of  gravity  the  damping  terms  evaluate  as: 

Hdampi  =-2C  *Bho*Vbar  1 *1  V bar  1 1 *Ldbar  1 * (L/2-X eg)  *Ws  1 
f 

(11-28) 

Hdaap2=-2C*Bho*Vbar2* IVbar2 | *Ldbar2* (L/2+Xcg) *Ws2 

(11-29) 


e.  Planing  Forces 

Initial  time  domain  validation  studies  conducted 
with  the  model  developed  from  the  equations  of  this  chapter, 
and  the  Loads  and  Motions  program,  disclosed  a serious 
discrepancy  between  the  two  models  in  natural  pitch 
frequency.  Specifically,  the  Loads  and  Motions  program 
predicted  a natural  pitch  frequency  of  3.2  rad/sec  at  10 
lets,  and  5.4  rad/sec  at  20  kts,  compared  with  a constant 
pitch  frequency  of  3; 3 rad/sec  for  the  linear  model.  This 
was  taken  to  indicate  the  presence  of  a velocity-dependent 
source  of  theta  feedback.  Examination  of  ref  disclosed 

that  Leo  and  Boncal  had  indeed  utilized  an  added  mass  effect 
which  had  the  effect  of  a planing  force,  in  their 
modifications  of  the  Loads  and  Motions  program. 

A body  in  motion  in  a fluid  experiences  a 
lifting  force  proportional  to  the  square  of  the  velocity  and 
directly  proportional  to  the  angle  of  attack,  or: 

2 

Hplan  =-C*Theta*Rho*V»s/2 

la 


where  S is  the  surface  area  of  the  body,  or: 


(11-30) 
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s 


= (HsUWs2)  *L/2 


(11-31) 

and  C is  the  lift-curve  slope.  This  term  is  not  directly 
la 

available,  and  is  a complex  function  of  the  sidewall 

geometry.  A value  of  1.5  gave  good  agreement  with  the 
natural  frequency  of  the  Loads  and  Motions  program,  but  a 
less  empirical  approach  to  this  value  is  recommended  for 
future  studies. 

The  planing  force  is  assumed  to  act  at  the  aft 
quarter-length  point. 

3 . Sum  of  Mpm ents 

The  forces  described  above  are  multiplied  by  their 
respective  moment  arms  as  shown  in  Fig  2 to  yield  the 
pitching  moments  about  the  center  of  gravity.  The  angular 
acceleration  about  the  C.G.  is  then: 

Alfa  = (Ppres+Pbuoy+Pseal+Pdamp+Pplan) /lyy 

(11-32) 

a.  Pressure  Moment 

The  pressure  exerts  a moment  about  tne  center  of 
gravity,  equal  to: 
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Ppres  =-Hpres*  (Xcp-Xcg) 


(II-33) 


b.  Buoyant  Moments 

The  buoyant  forces  due  to  the  immersed  vclumes 
fore  and  aft  of  the  center  of  gravity  exert  a pitching 
mcment  along  their  moment  arms,  which  are  approximately  (for 
small  angles)  the  centroid  of  the  rectangle  from  the  center 
of  gravity  to  the  bow  and  stern: 

Pbuoyl  =-. 5*  (L/2- Xcg) ♦Hbuoyl 

(11-34) 

Pbuoy2  =. 5 * (L/2+Xcg)  *Hbuoy2 


(11-35) 


c.  Bow  and  Sternseal  Moments 

The  bowseal  acts  along  a moment  arm  that  is 
formed  by  the  distance  from  the  C-G.  to  the  centroid  of  tiie 
wetted  area  Aseal,  or: 

Pseall  =-Hseal1* (L/2-Xcg-Xseal 1/2) 

(11-36) 

The  sternseal  force  acts  on  a moment  arm  that  is 
formed  by  the  distance  from  the  C.G.  to  the  centroid  of  tae 
sternseal  wetted  area,  or: 
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Fseal2 


= Hseal2*  (L/2+Xcg+Xseal/2) 


(11-37) 


d.  Damping  I'lcments 


The  damping  moments  in  pitch  are  the  damping 
forces  in  heave  multiplied  by  their  respective  moment  arms, 
two-thirds  the  distance  from  the  C.G.  to  the  bow  and  stern 
respectively : 

Pdampi  =- (2/3) *(L/2-Xcg)  ♦Kdampi 

(11-38) 


Pdamp2  = (2/3)  » (L/2  + Xcg)  *Kdamp2 


(11-39) 


e.  Planing  Moments 


The  planing  moment  is  the  planing  force  acting 
at  the  aft  quarter-length  point,  or: 

Pplan  = (L/4+Xcg)  »Hplan 


(11-40) 


2.  INITIAL  CONDITIONS 


As  in  the  Loads  and  Motion  program,  initial  conditions 
must  be  carefully  selected  to  ensure  steady  state  conditions 
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at  the  beginning  of  a run.  Since  speed  does  not  enter  into 
this  program,  one  of  the  major  problems  of  the  Loads  and 
Motions  program  is  obviated.  However,  the  craft  in  this 
model  is  particularly  sensitive  to  pitch,  and  once  excited 
in  this  mode,  is  very  lightly  damped. 

Selection  of  initial  conditions  proceeds  as  follows: 


1.  pitch  rate  and  vertical  velocity  are  defined  as  zero; 

2.  a value  of  Pfabar  is  selected  such  that  plenum  air  mass 
flow  rate  is  set  to  zero; 

3.  a value  of  draft  is  calculated  such  that  weight  and 
pressure  forces  are  balanced  by  buoyant  and  planing 
forces ; 

4.  a value  of  Theta  is  calculated  such  that  pressure 
moments  are  balanced  by  buoyant  and  planing  moments; 

5.  given  the  Plenum  volume  defined  by  the  equilibrium 

values  of  draft  and  Theta,  and  the  equilibrium  value  of 
Pbbar,  the  Plenum  air  mass  may  be  calculated 

adiatatically . 


More  detail  on  the  exact  procedures  for  the  selection  of 
equilibrium  conditions  may  be  found  in  subsequent  chapters 
dealing  with  the  computer  simulation  of  this  model. 
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Figure  3 - DETAIL  OF  BOW  AND  STERK3EAL 


6 


III. 


LINEAR  HEAVE-PITCH  MODEL 


The  non-linear  equations  of  motion  developed  in  the 
preceding  chapter  form  the  basis  of  a linearized  heave-pitch 
model,  valid  for  disturbances  about  its  equilibrium  point, 
defined  as  that  point  aoout  which  all  rates  are  z ero . The 
linear  model  permits  the  use  of  a variety  of  powerful 
techniques  for  parametric  analysis,  such  as  root  locus;  it 
furthermore  simplifies  the  design  of  an  automatic  control 
system,  which  could  serve  to  further  limit  excursions  about 
the  operating  point. 

A.  TAYLOR'S  SERIES  EXPANSION 

For  small  excursions  about  an  operating  point, the 

non-linear  relationships: 

f =F(x,y) 

g =G(x,y) 

are  approximately: 

f +df  =F  (X , V)  + ( dF/d  X)  *Dxt  (d?/d  V)  *Dv 

0 0 0 

g+  dg  =G(x  ,v  ) t (dO/dx) *Dx+ (dG/dy ) *Dy 

0 0 ' 0 

where  the  derivatives  of  F and  G with  respect  to  x and  y are 

evaluated  at  the  operating  ooint,  x and  y.  Cancellinq  the 

0 0 

equilibrium  point  from  both  sides  and  dividing  through  by 
time  yields; 
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d/  it)  =(jr/ix)*  1/  (X)  * (d//dy)  ^ d/  (/) 

/dt  /at  /dt 


d/  (q) 
/dt 


= (d^/.ix)*  i/  i'x)  { L,/djj  * d/  (i) 
/at  /Jt 


Appiyiiiq  the^e  pr  i ncijljij  tJ  ca-i  non- o i tn.; 
pr <=  cea iu  ] chaptei',  thv^  ot^iiy  state  opei  atin;  pciat  is 
defined  as  dcfore  tor  th>-j  eq  ai  inr  i am  vdiues  of  tiia  stare 
/ jriabies . 


d.  PLdNUK  JaUGh  JZHI.’AriVGG 

Plenjii  qaaq=  uressuLe  has  o.aei.  dedi:.  ‘i  as; 

PLtLar  =?.)-Pd 


(ZI---)) 


and  since  ?a  is  a coastdjit: 
di'tcar  =d?o 

Pn,  the  pleauui  totsj.  -ressure,  aas  eotn  Itfinc..!  as: 


Po  (0) 


= ? i*  (.id  (d)  / (Zuoa  ‘ / 3 (J)  I 


bd  II  n i 


(ri-d) 

Thecefore  eXrdnaiaq  around  th>-.  eq  ui  1 Lo  ri  a.r.  vafios  of  ^s(O) 
and  d b ( 0 ) : 


t b (0)  ♦ d r !)  = P \ * ( d b ( J ) / ( R hod  * V b (d)  ) 


.d  ;i  n i 


4 r 


+ (Gamma/Mb  (0)  ) *Pa*  (Mb  (0)  / (fihoa^Vb  (0)  ) 


Gamma 

- (Gamna/Vb (0) ) *Pa* (Mb (0) / (Rhoa*Vb ( 0) ) *dVb 

(III-1) 

or,  after  substituting  the  definition  of  Pb(0)  and 
cancelling  the  eguil^-brium  values  from  both  sides: 

dPb  =Gamma*  (Pb  (0) /Mb  (0)  ) *dMb 

-Gamma*  (Pb  (0)  /Vb  (0)  ) *dVb 


(III-2) 

Vt (0)  may  be  further  expanded, since  it  has  been  defined  in 
the  preceding  chapter  as  the  empty  plenum  volume  less  the 
net  submerged  volume,'  following  the  definition  of  ref  W: 

Vb(0)  =Vn-Ab*  (z(0)  + Zs+Xcg*Theta  (0)  ) 

(II-7) 

and 

dVb  =-Ab* dz-Ab*Xcg*dTheta 

(III-3) 

Substituting  the  expansion  of  incremental  volume,  Fguation 
(111-3) , into  the  expansion  for  incremental  pressure. 
Equation  (III-2) , and  dividing  by  differential  time  yields: 

d/  (Pbbar)=Gamma*  (Pb  (0) /Mb(0)  ) * d/  (Mb) 

/dt  /dt 
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♦ d u d * A L * ( ? b ( J ) / V i;  ( 3 ) ) *■ 


{-) 

■/  (Ihetd) 

/It 

( Ill-d) 

Fct  ;;iOL^  coinpdct  nordtioi.  t'ue  foiljwiu.j  iniiuencd 
coc  f tijidntii  are  .sacn  that: 

L J D Z d ,ii  I!1  d ^ A L •*“  ( P u ( 3 / / ’.'  d ( .3  3 ) ; 


J/ 
/at 

+ d a n ui  d-i"  A b ^ X c g * I p u ( J ) / '/  b ( 0 ) ) ’* 


( li  I-T) 


l)  ab  t ii=Gd  in.Jd  » Ad  »Xc  ( P o i 3 ) / V b {')  ) ; 


( 


■J) 


3 a b :n  d = jd  3 1 d =*  ( ? d ( 3 ) / -i  d ( 0 ) ) ; 


and  th-j  tresdure  1 1 f i dL'aa t ia  1 -,-^uitidii  d-'co.n.'-d : 


d/  (Prbar)=Dpbz*  d/  (z) 
/at  /.it 


^•Dddth*  li/  (Thriti) 
/.it 


♦ OdJ.nt^  1/  f/!D) 

/it 


(Tir-7) 


llir-.) 


t'LzNJ.t  A 13  .1A.S3  FlC/.' 


d A r d D d I V A i'  I V d G 


F 


The  plenum  air  mass  flow  rate  has  been  defined  in  the 
preceding  chapter  as; 


d/  =Rhoa*  (Qin  (0) -Qout  (0)  ) 

/dt 


and 

d(  d/  (Mb)  )=ahca*dQin-Rhoa^dQout 
/dt 

Qin  has  been  defined  as: 

Qin  (0)  =N*  (QiO-Pbbar  (0) ) 

so  that 

dQin  =-N*dP.bbar 


(11-12) 


(III-9) 


(11-10) 


(III-10) 


and  Qout  has  been  defined  as; 

1/2 

Cout(O)  =Cn*Al*  (2*Pbbar  (0) /Rhoa) 

(11-11) 

and 


-1/2 

dQout  = (Cn*Al/Rhoa)  * (2*Pbbar  (0)  /Rhoa)  *dPbbar 

(III-1  1) 

Combining  the  two  expressions  for  incremental  volumetric 
flew  into  and  out  of  the  plenum.  Equations  (111-10,11),  with 


51 


the  expression  for  ior  incremental  mass  flow  rate.  Equation 
(III-9) , and  dividing  by  differential  time  yields: 

2 -1/2 
d/  2 (Mb)  =- (N*fihoa+Cn*Al*  (2*Pbbar  (0) /Ehoa)  ) ♦ d/  (Pbtar) 
/dt  /dt 

(III-12) 

The  influence  coefficients  may  be  defined  by  substituting 
Equation  (III-8)  into  Eg uat ion (III-1 2)  and  equating  terms: 

-1/2 

Dmbz  =- (N*Rhoa+Cn^Al*  (2*Pbbar (0) /Rhoa) ) *dpbz 

(III-13) 


-1/2 

Dmbth  =- ( N*ahoa+Cn*Al* (2*Pbbar (0)  /Rhoa) ) * dpbt  b 

(III-14) 


-1/2 

Dmbmb  =- (N*ahoa+Cn*Al* (2*Pbbar (0) /Rhoa) ) *dpbmt 


(III-15) 


Thus  the  linearized  differential  equation  for  plenum  mass 
flow  rate  may  be  written: 


2 

d/  2 (Mb)  =Dmbz*  d/  (z) 
/dt  /dt 


♦Dmbth* 

(Theta) 

/(it 

♦ Dmbmb* 

(Mb) 

/dt 

(III-16) 
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D.  HEAVE  DEfilVATIVES 


The  heave  equation  has  been  defined  as; 

2 

d/  2 (z)  = (Hpres+Hbuoy  + Hseal+Hdamp  + Hplan+H) /a 

/dt 

(11-13) 

Danifing  forces  are  neglected  in  the  linear  model. We  shall 
consider  the  various  components  of  heave  force  separately, 
and  superimpose  the  components  upon  completion  of  all 
derivations. 

1 . Heave  Fores  ^ue  To  Pressure 

The  heave  force  due  to  pressure  has  been  defined  as; 
Hpres  =-Ab*P.nbar 


(11-14) 

By  inspection,  its  ioifluence  coeff icienrs  become; 

Dhpz=-Ab*Dpbz 


Dhpth=-Ab*Dpbth 


1 Dhpmb=-Ab*Dpbmb 

i. 

I 


1 

1 


(III-17) 


(III-18) 


(III-1  9) 
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2 • H€avg  Forces  Due  To  Buoyancy 

The  buoyant  heave  force  equation  from  the  preceding 
chapter  has  been  defined  as: 

Hbuoy  =Hbuoy 1 +Hbuoy 2 

=-2*Rho*G*Ws  1*  (L/2-X.cg)  * (Ld-  (L/2-Xcg)  *Theta/2) 

-2*Rho*G*Ms2*  (L/2  + Xcg)  ♦ {Ld+  (L/2+Xcg)  «Tneta/2) 

(11-15,16,  17,13) 

which  may  be  rearranged: 

Ebuoy= 

-2*Rho*G*<Wsl*  (L/2-Xcg)  + ws2*  (L/2+Xcg)  ) * (z  (0)+Zs) 

2 2 

+ 2*Rho*G*  (Ms  1*  (L/2-X.cg)-Ws2*  (L/2+Xcg)  ) *Theta/2 

and  the  influence  coefficients  may  be  written  by  inspection: 
Dhbz  =-2*Rho*G* (Ws 1* (L/2-Xcg) + Ws2*  (L/2+X eg) ) 

(III-20) 

2 2 

Ehbth  =Rho*G*  (Wsl*  (L/2-Xcg)  -Ws2*  (L/2  + Xcg)  ) 

(III-21) 


3 . He ayg  Fo£c « due  ^ Bow  and  S ternseal 


■i 

'1 


1 


1 

) 

i 


i 


I 
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The  expression  for  the  bowseal  heave  force  has  been 
defined  as: 

Hseall  (0}  =-Midth*Xseal  1 (0)  *Pbbar  (0) 

(11-21) 

which  for  saall  perturbations  in  Xseall  and  Pbbar  becoaes 
dHseall  =-Hidth*Pbbar  (0) »dXseal1 

-Hidth*Xseal1  (0)  *dPnbar 


[111-22) 

The  expression  for  Xseall  has  been  defined: 

Xseall  (C)  =1.65*  (z(0) +Zs- (L/2-Xcg) *Theta) 

(11-20) 

which  for  saall  perturbations  around  z(0)  and  Tneta  (0)  is 
approximately 

dXseall  =1.65*dz-1.65* (L/2- Xcg) *dTheta 

(III-23) 

Substituting  these  relations  along  with  the  relations  for 
small  perturbations  in  Pbbar  due  to  z, Theta  and  Mb  from 
Equations  (111-5,6,7)  yields; 

dHseal  1 =-Midth*(1.65*pbbar(0)  +Xseal  1 (0)  *Dpbz)  *dz 

-Width* (- 1 .65* (L/2-Xcg) *Pbbar  (0/ 

♦Xseall (0) *Dpbth) *dTheta 
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- Hidth*Xseai  1 (0)  *Dpbinb*dHb 


(III-24) 

Th€  influence  coefficients  then  become: 

Dhszl  =-Width*(1 . 6 5*?bbar  (0)  +Xseall  (0)  ♦Dpbz) 

(III-25) 


Dhsthi  =-Widt  h*  (-1.6  5*  (L/2-Xcg)  *Pbbar  (0) 

+Xseal1  (0)  *DFbth) 


(111-26) 


Dhsmb  =-Hidth*Xseal1  (0)  *Dpbmb 


(III-27) 


The  expressipn  for  the  stern  seal  has  been  defined: 
Hseal2  =-3ealp*Wid-ch*Xseal2  (0) 

(11-23) 

Xseal2(0)  =1.65*  (Ld  (0)  + (L/2  + Xcg)  *tan  (Theta  (0)  ) ) 


(11-22) 

Thus  for  small  perturbations  in  draft  and  Theta: 
dHseall  =- 1 .65*sealp*Midth*dz 

-1.65*  (L/2  + Xcg)  *Sealp*Width*dTheta 

(III-26) 

Thus  the  sternsaal  influence  coefficients  become: 

i 

I 


5b 


=-1 .65 ♦Seal?*Width 


(III-29) 


Cbsth2 


=-1  .65*  (L/2  + Xcg)  *Sealp*Width 


(III-30) 


to  Planing 


The  planing  force  is  linear  in  Theta  , or: 


Hplan 


=-C  *fiho*7*  (Ws  1+Ws2)  *L*Theta 
la 


(11-30,31) 

and  the  influence  coefficient  may  be  written  by  inspection: 


Dhplth 


=-C  ^Rho*V*  (Ms  UWs2)  *L 
la 


(III-31) 


5 . Sum  Cf  Forces 


The  forces  previously  derived  may  be  superimposed  so 


that ; 


2 2 
d/  2 (z  (0)  ) +d  ( d/  2 (z)  ) = 
/dt  /dt 
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(Hpres  (0)  ♦Hbuoy  (0)  ♦Hs  eal  (0)-mpian(0))/a 
* (Dhpz+Dhbz  + Dhs 2)  *dz/M 
(Dhpth+Dhbth+Dhsth  + Dhplth)  *dTheta/M 
(Dhpfflb+Dhsmb)  *dMb/M 

(III-32) 

The  influence  coefficients  aay  now  be  defined: 

Dzz=  (Dhpz+Dhbz+Dhsz)  /fl 


(III-33) 

Dzth=  (Dhpth  + Dhbth+Dhsth+Dhplth)  /ii 


Dzmb=  (Dhpab  + Dhsmb) /M 


(111-3*4) 


(III-35) 


Heaoving  the  equilibrium  values  and  dividing  through  by 
differential  time  yield  the  differential  equation  of  motion 
in  the  z-direction: 


3 

d/  3 (2)  =Dzz^  d/  (z) 
/dt  /dt 


♦ Dztn*  d/  (Theta) 

/dt 


♦ Dzmb*  d/  (Mb) 
/dt 


(III-36) 


2.  PITCHING  HCMENT  DERIVATIVES 
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The  pitching  aoaent  equation  is 


d/  2 (Theta)  = (Ppres+Pbuoy  ♦Pseal+Pdaap+Pplan) /ly y 
/dt 

(11-32) 

As  in  the  previous  case  of  heave  forces,  damping  moments  are 
neglected  in  the  linear  model;  the  various  moments  are 
considered  separately  and  superimposed  after  derivation. 

• Pitch  Moment  due  to  Pressure 

The  pressure  moment  has  been  defined  as 
Ppres  = (Xcg-Xcp)  *Hpres 


(11-33) 

The  influence  coefficients  may  be  written  by  inspection: 

Dppz=  (Xcg-  Xcp)  *Dhpz 


Dppth=  (Xcg-Xcp)  *Dhpth 


Oppmb=  (Xcg- Xcg)  *Dhpmb 


(III-37) 


(III-38) 


(III-39) 


2.  Buoyant  Moments 


The  huoyant  moment  has  oeen  defined: 
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r 


I 

f 

! 


Pbuoy  =Pbuoy  1 ■♦•Pbuoy2 

=-  (L/2-Xcg)  ^-Hbuoy  1/2+  (L/2  + Xcg)  »Hbuoy2/2 

2 

= Rho*g*  (Wsl*  (L/2-Xcg)  * (Ld-  (L/2-Xcg)  *2hata/2) 

2 

-Ws2*  (L/2  + Xcg)  * (Ld+  (Ld+  (L/2  + Xcg)  =»<Theta/2)  ) ) 

(11-34,35) 

and  the  influence  coefficients  aay  be  defined: 

2 2 

Dpbyz  =Rho^G*  (Msl*  (L/2-Xcg) -Ws2* (L/2  + Xcg)  ) 

(III-40) 

3 3 . 

Dpbyth  =-Rho.«G*(Ws1*  (L/2-Xcg)  +Ws2*  (L/2  + Xcg)  ) /2 

(III-41) 

3 • Sow  and  5 ternseal  Moments 

The  fowseal  contributes  a moment 
Pseall  =-Hseal  (L/2-Xseal 1/2-Xcg) 

(II-3o) 

which  for  small  disturbances  , is  approximately 
dPseall  =- (L/2-Xseal 1 (0)/2-Xcg) ♦dHseall 
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J 


♦Hseall (0) *dxseal 1/2 


(III-42) 

Substituting  the  definition  cf  dXseall,  Equation  (III-23) , 
and  the  influence  coefficients  of  Hseall,  Equations 
(111-25,26, 2'/)  yields: 

dEseall  =(-  (L/2-Xseal  1 (0)  /2-Xcg)  *Dhsz1 

♦ 1.65*Hseal1  (0)  /2)  *dz 
♦ (-  (L/2-Xseal1  (0)  /2-Xcg)  ♦Dhsthi 
- 1.65  (L/2-Xcg)  *Hseal1  (0)  /2)  ^dTheta 
- (L/2-Xseal1  (0)  /2-Xcg)  *Dhsint*dMb 


The  influence  coefficients  may  now  be  defined: 


Dpszl 


=-  (L/2-Xseal1  (0)/2-Xcg)  ^Dhszl 
+ 1.65*Xseal  1 (0)  /2 


(III-43) 


Dpsthi 


= - (L/2-Xseal1  (0) /2-Xcg)  ♦Dhsthl 
- 1.65*(L/2-Xcg)  *Hseal  1 (0)  /2 


Cnsmb 


=-  (L/2-Xseal  1 (0)  /2-Xcg)  ♦Dhsab 


(III-44) 


(III-43) 


(III-46) 


The  sternseal  contributes  a moment: 
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Fseal2  (0) 


= Hseal2*(L/2  + Xseal2  (0)  /2+Xcg) 


(11-37) 

which  for  small  disturbances  is  approximately: 
dPseal2  = (L/2+Xseal2  (0)  /2  + Xcg)  v-iHseal2 

+Hseal2 (0) *dXseal2/2 

(III-47) 

and  substitution  from  the  sternseal  force  influence 
coefficients.  Equations  (111-29,30),  yields; 

dPseal2  = ( (L/2  + Xseal2  (0) /2  + Xcg) *Dhsz2 

+ 1.65*Hseal2  (0)  /2)  *dz 

+( (L/2+Xseal2 (0) /2+Xcg)  *Dhsth2 

+ 1.65*(L/2+Xcg) *Hseal2 (0) /2) +dTheta 

(III-43) 

and  the  stern  seal  pitch  influence  coefficients  may  be 
defined: 

Dpsz2= (L/2+Xseal2(0)  /2  + Xcg)+Dhsz2 
+ 1.65*Hseal2  (0)  /2 

(III-49) 

Dpsth2  = (L/2  + Xseal2  (0) /2  + Xcg) *Dhsta2 
+ 1.  65*  (L/2+Xcg)  *Kseal2  (0)  /2 

(III-50) 

• Pitch  Mgment  due  to  Planing 
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The  flanicg  moment  is  linear  in  Theta; 


Pplan  = (L/4+Xcg)  *Hplan 

and  the  influence  coefficient  is; 
Dpplth  = (L/4+Xcg) *Dhplth 


(11-40) 


(III-51) 


5 . Sum  foments 


The  individ-ual  moment  contr  ibutions  may  now  be 

summed ; 

Alfa(O)  ♦dAifa  = (Ppr6s(0)>Pbuoy  (0)  ♦Pseal  (0)  +Pplan  (0)  ) /ly y 

♦ (Dppz+Dpbyz  + Dpsz)  *dz/Iyy 
♦ (Dppth+Dp byth+Dpst h+Dpplth) *DTheta/Iyy 
♦ (Dppmb+Dpsmb)  *dMb/Iyy 


The  influence  coefficients  may  now  he  defined 

Dpz.=  (Dppz  + Dpby  z+Dpsmb) /I  yy 


(III-52) 


(III-53) 

Dpth=  (Dppth+Dpbyth+Dpsth  + Dppith) /lyy 


(III-54) 
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Dpab=  (Dppob+Dpsmb)/Iyy 


(III-55) 


and  the  differential  equations  of  motion  for  angular  rate  of 
acceleration  nay  be  written,  after  dividing  through  by 
differential  time  and  cancelling  the  equilibrium  values; 


3 

d/  3 (Theta)=Dpz*  d/  (z) 
/dt  /dt 


+Dpth* 

/dt 

(Theta) 

•fDpmb* 

/dt 

(Mb) 

(III-56) 


F.  STATE  SPACE  EEPRSSENT ATION 

Fig  4 is  the  linear  model  as  derived  in  this  chapter 
presented  in  matrix  form  as; 

X =A*x 

Examination  of  the  eigenvalues  of  the  9-by-9  matrix  shows 
four  eigenvalues  at  the  origin  (0±0j)  which  represent  the 
integrated  output  equations  tor  z.  Theta,  Mb,  and  Phbar;  the 
remaining  eigenvalues  consist  of  a complex  pair,  located  in 
the  vicinity  of  *30±30j,  and  a real  eigenvalue  in  the 
vicinity  of  -.5,  as  in  ref  W.  and  an  additional  complex 
pair  near  0±4j,  representing  the  very  lightly-damped  pitch 
response . 
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IV.  TIHE  DOMAIN  VALIDATION 


The  linear  heave-pitch  XR-3  model  was  validated  against 
the  modified  Oceanics  XR-3  Loads  and  Motions  Program  cf  ref 
prior  to  investigations  in  the  frequency  domain.  The 
linearized  program,  which  appears  at  the  end  of  the  text, 
was  simulated  using  the  Continuous  System  Modeling 

Program,  version  III'.  This  is  basically  a FORTRAN  language 
with  extremely  powerful  variable-step  integrators  and 

versatile  output  structure.  Readers  familiar  with  FORTRAN 
should  find  little  difficulty  following  the  CSMP  III 

programs  developed  for  this  simulation;  for  details  of  this 
simulation  language,  see  ref  M- 

The  Time  Domain  Program  consists  of  three  sections: 

1.  Initial  Conditions;  this  subsection  identifies  the 
equilibrium  operating  point  for  a given  set  of  craft 
physical  parameters. 

2.  Linearization;  the  partial  derivatives  of  the 

non-linear  state  equations  derived  in  the  preceding 
chapter  are  evaluated  at  this  particular  operating 
point. 

3.  Dynamics;  the  weight  and  C.G.  acceleration  are 

perturbed  at  time  zero  by  a step  weight  change,  and  the 
resultant  linear  state  trajectories  in  time  are 

computed.  Note  that  the  linear  model  is  not 

self-starting;  an  initial  C.G.  acceleration  must  also 
be  applied. 
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L 


A 


INITIAL  CONDITIONS 


Initial  conditions  involve  the  solution  of  the  | 

non-linear  system  equations  of  Chapter  II  for  the  j 

equilibrium  operating  point  defined  at  the  end  of  the  ; 

chapter. 

1 Presgare 

t 

A value  of  plenum  pressure  is  selected  so  that 
plenum  mass  flow  rate  is  set  to  zero,  or,  from 
Equations(II-10,11) : 

Cin  =Qout 

V2 

j giO-Pbbat  = (Cn*Al/N)  * (2*Pbbar/Rhoa) 

i 

■ (IV- 1) 

I 

which  is  a quadratic  in  Pbbar: 

2 2 2 
Pbbar- (2*QiC+ (Cn*Al/N)/Rhoa) *Pbbar+Qi0=0 

(IV-2) 

The  larger  solution  of  this  quadratic,  correspondi ng  to  the 
physically  unrealizable  condition  of  reverse  flow  through 
; the  fans,  is  rejected. 

■ 2.  £raft  and  git^h  Angle 

I 

I 

The  solution  for  equilibrium  values  of  draft  and 
pitch,  given  the  previously-computed  value  of  Pbbar, 
requires  an  iterative  technique,  due  to  their  close 


n 


coupling.  An  initial  value  of  draft  is  first  computed 
assuming  zero  pitch  and  no  bow  or  stern  seal  contribution. 
The  ncn-linear  state  equations  are  then  solved  for  residual 
forces  and  moments  due  to  all  sources.  A very  small 
perturbation  in  draft  and  pitch  angle  of  one  percent  with  a 
minimum  default  magnitude  of  .001  is  applied  and  tne 
differential  change  in  residual  forces  and  moments  with 
respect  to  draft  and  Theta  are  computed.  Eairstow's  method 
for  the  zeroes  of  functions  or  two  variables,  as  described 
in  ref  , is  used  to  drive  the  residuals  to  zero.  This 
particular  problem  is  somewhat  slow  to  converge;  convergence 
is  enhanced  by  amplifying  the  computed  changes  in  draft  and 
pitch  by  1.9,  but  after  two  hundred  iterations  residuals  on 
the  order  of  ten  Ihf  and  ten  ft-lbf  may  still  be  observed- 
These  residual  forces  and  moments  are  applied  to  the  initial 
values  of  22dot  and  Alfa. 

3.  Plenum  Air  Mass 


The  equilibrium  values  of  draft  and  pitcn  angle 
define  a unique  plenum  volume; 

Vb  =Vn-Ab*  (Ld+Xcg*tan  (Theta)) 


(II-7) 

which,  when  combined  with  the  equilibium  pressure,  yields  a 
unique  adiatatic  plenum  air  mass,  from  Equation  (II-8J  : 

1/Gamma 

MbO  =Rhoa!»‘ Vb*  (Pbbar/Pa) 


(IV-3) 


The  solutions  to 
equilibrium  values  of  Pb (0)  , 


these  equations  define  the 
Ld(0),  Theta(O),  and  Mb(0) 
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about  which  the  partial  differential  equations  of  the 
preceding  chapter  may  be  evaluated. 

B.  LOADS  AND  MOTIONS  PROGRAM  OPERATING  CONDITIONS 

The  XR-3  Loads  and  Motion  program  was  run  at  6500  lb 
weight,  C.G.  S.825  ft.  ahead  of  the  aft  transom  at  20  kts. 
Simulation  time  was  30.0  sec.  for  the  purpose  of  obtaining 
reasonablj  quiescent  values  of  pitch  angle,  draft,  and 
Plenum  pressure.  In  order  to  obtain  identical  operating 
conditions,  it  was  necessary  to  adjust  various  constants  in 
the  linear  model  until  agreement  was  reached  with  the 
steady-state  operating  conditions  of  the  Loads  and  Motion 
program.  Specifically,  these  parameters  were; 

1.  seal  leakage  area,  Al,  adjusted  until  the  two  programs 
agreed  in  Plenum  pressure 

2.  Planing  lift-curve  slope  coefficient,  C , adjusted 

la 

until  the  two  programs  agreed  in  natural  pitch 


frequency 

9 

3. 

minimum 

sidewall  width,  WsO, 

adjusted 

until 

t he 

two 

programs 

agreed  in  draft; 

4 . 

Residual 

pitching  moment  ,Mf0, 

adjusted 

until 

t he 

two 

programs 

agreed  in  steady-state 

pitch  angle. 

The  resul 

ts  of  these  trial-and-e 

rror  iterations 

yielded 

the 

following 

table  pf  values: 
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TABLE  OF  OPEBATJNG  PARAHETEES 


;,JNEAR  SYSTEM 

LOADS  AMD  MOTIONS 

AfFECTED 

STEADY-STATE 

STEADY-STATS 

PARAMETER 

VfEIABLE 

VALUE 

VALUE 

2 

2 

Al=.390 

Ptbar 

24.85  Ibf/ft 

24.84  Ibf/ft 

C=1.5 

5.02  rad/sec 

5.46  rad/sec 

la 

F 

Ws0=.27 

Id 

8.497  in. 

8.55  in. 

Mf0=-1875 

Theta 

1.21  deg 

1.199  deg 

C.  COMPASISCN  OF  RESULTS 

The  Loads  and  Motions  program  was  run  3.0  secs  with  tne 
initial  conditions  as  obtained  above,  but  with  a weight  of 
6400  lbs.  The  linear  model  was  simultaneously  run  with  a 
100  lb  step  weight  reduction;  data  from  the  Loads  and 
Motions  program  was  subseguently  entered  into  the  CSMP  III 
linear  program  in  tabular  form,  for  simultaneous  plotting 
and  error  computation.  For  the  purposes  of  comparison,  error 
was  defined  as  the  difference  between  the  data  from  the  two 
programs,  divided  by  the  average  value  of  the  two.  This 
procedure  significantly  reduced  the  probability  of 
inadvertent  division  by  zero,  as  compared  with  normalizing 
the  difference  with  respect  to  one  system  or  the  other. 

^ • Pi§S!i5  Pressure 

Plenum  pressure  and  Plenum  pressure  error  (percent) 
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plotted  vs.  time  ippear  in  Figs  5 and  6.  The  two  curves  are 
in  quite  good  general  agreement,  with  the  magnitude  of  error 
being  on  the  order  of  several  parts  per  thousand.  The 
linear  model  exhibits  a greater  pressure  reduction  in  the 
initial  pressure  transient,  than  does  the  Loads  and  -lotion 
results,  and  this  transient  is  more  heavily  damped  in  the 
linear  model.  One  explanation  for  this  difference  would 
appear  to  be  the  variable  leakage  area  of  the  Loads  and 
Motion  program.  Pitch  appears  to  couple  more  strongly  into 
pressu're  in  the  linear  model;  however,  the  magnitude  of 
pitch  is  significantly  great-ar  in  the  linear  model  due  to 
the  absence  of  significant  damping  terms.  Finally,  the 
pressure  in  the  Loads  and  Motion  program  is  tending  toward  a 
slightly  higher  value,  further  supporting  the  supposition 
that  leakage  area  has  decreased  in  the  Loads  and  motion 
program. 


C.G.  Acceleration 


Figs  7 and  8 show  C.G.  acceleration  and  C. G. 
acceleration  error  vs.  time.  They  generally  exhibit  similar 
characteristics,  but  due  to  the  larger  initial  pressure 
reduction  of  the  linear  model,  there  is  a more  pronounced 
overshoot  in  C.G.  acceleration  in  the  linear  model.  In 
addition,  the  stronger  pi tch -pressur e coupling  in  the  linear 
model  also  translates  as  stronger  pitch-C.G.  acceleration 
coupling.  The  errors  appear  as  relatively  large 
percentages;  however,  these  are  the  differences  of  two  very 
small  numbers,  and  should  be  regarded  cautiously. 


3.  Draft 


10. 


Draft  and  draft  error  vs.  time  appear 
As  has  been  previously  verifed  by  the 


in  Figs  9 an-l 
Ger ba/Ihaler 
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model,  draft  is  very  well  approximated  by  a linear  model. 
For  this  particular  disturbance,  it  appears  to  be  slightly 
less  damped  and  tending  toward  a lower  steady-state  value  in 
the  linear  mcdel.  Little  evidence  of  pitch  coupling  is 
apparent,  and  errors  are  on  the  order  of  parts  per  thousand. 

4 . Pitch  Rate 

Pitch  rare  and  pitch  rate  error  are  shown  plotted 
vs.  time  in  Figs  11  and  12.  The  linear  model  exhibits  a 
pitch  rate  of  approximately  twice  the  magnitude  of  the  Loads 
and  Motion  program;  this  may  ostensibly  be  attriouted  to  the 
sharper  pressure  transient  previously  observed;  with  tne 
center  of  pressure  fprward  of  the  C.G.  in  both  models,  this 
would  induce  a higher  initial  angular  acceleration  in  tne 
linear  model.  In  addition,  the  linear  model  is 
significantly  less  damped  in  pitch;  therefore  this  larger 
magnitude  angular  acceleration  would  build  to  a larger 
angular  velocity,  and  persist  for  a longer  period  of  time  in 
the  linear  mcdel.  The  frequencies  compare  favorably,  5.0U 
rad/sec  for  the  linear  model  and  5.28  for  the  Loads  and 
Mction  program. 

5 • Pitch  Angle 

Figs  13  and  14  display  pitch  angle  in  degrees  and 
pitch  angle  error  vs.  time.  Both  show  basically  identical 
characteristics;  i.e. , an  initial  pitch  down  induced  by  the 
pressure  drop,  and  an  exponentiai  decay  to  a lower  pitch 
angle  induced  by  the  draft  change.  The  magnitude  of  pitch 
angle  change  reflects  the  larger  angular  velocity  of  the 
linear  mcdel,  and  the  difference  between  the  damping  of  the 
Loads  and  Motion  program  and  the  linear  model  is  quite 
dramatic.  The  magnitude  of  error,  however,  is  on  the  order 
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Figure  14  - PITCH  ERROR  VS.  TIME 
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V.  DOMAIN  DEVELOPMENT 


This  chapter  develops  the  techniques  for  frequency 
response  of  a general  linear  system  by  complex  matrix 
inversion.  The  ocean  wave  equations  utilized  by  Booth  are 
reviewed  and  summarized,  and  the  forcing  functions  of  wave 
action  on  the  linearized  XR-3  model  are  developed  in  pitch 
and  heave.  Portions  of  Booth's  work  are  repeated  utilizing 
the  linearized  model,  and  the  results  contrasted  and 
compared , 


A.  LINEAR  SYSTEM  EREQOENCY  RESPONSE 


The  development  of  a linear  system  permits  the 
utilization  of  various  techniques  for  the  development  of 
frequency  response  characteristics  which  are  vastly  superior 
in  CPU  time  in  simulation,  than  the  methods  to  which  one  is 
limited  in  non-linear  models.  The  definition  of  the  state 
matrix  of  the  previous  chapter  defines  the  state  transfer 
matrix: 

-1 

H(s)  = I sI-A  I (s) 

(V-1) 

and  if  the  system  is  forced  by  a sinusoid,  s tends  toward  1w 
as  the  transients  decay  and: 


B4 


(V-3) 

9(w)  =tan I imag  (H  ( j w) /real  (H  ( jw)  j 

(V-4) 

where  G(w)  is  the  gain,  in  decibels,  with  respect  to  some 
arbitrary  reference  H(0)  , and  9 (w)  is  the  phase  angle  in 

the  complex  plane.  Thus  the  frequency  response  of  a linear 
system  tc  a given  sinusoidal  forcing  function  may  ne 
obtained  by  complex  matrix  inversion.  While  this  is  not  in 
itself  a particularly  economical  use  of  CPU  time,  it  is 
still  vastly  superior  to  a Fourier  transform.  This  is 
particularly  true  j.f  the  range  of  frequencies  is  wide  and 
variable-step  integrators  driven  to  small  step  sizes  which 
are  not  required  for  most  of  the  array  integration. 
Finally,  it  is  not  necessary  to  allow  the  system  to  attain 
steady-state  AC  conditions  pr  ar  to  computing  the  frequency 
response.  This  is  particularly  important  in  any  model  of  the 
XH-3,  which  exhibits  very  light  pitch  damping. 

B.  WAVE  CHAfiACTEHISTICS 

Booth  in  ref  provides  an  excellent  discussion  on  the 

characteristics  of  gravity-free  surface  waves  with  straight, 
infinitely  long,  parallel,  equally-spaced  wave  crests,  and 


I 
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constant  wave  height.  Such  a wave  advances  with  a velocity 
C which  is  often  referred  to  as  celerity,  in  order  to 
emphasize  the  fact  tiat  it  is  the  wave  front  which  advances, 
not  the  water  particles. 

The  regular  wave  may  thus  be  modeled: 

C =G/w 


(V-5) 


Lambda 


=2*Pi*C/w 


=2#Pi»G/w 

i 


(V-6) 

where  w is  the  incident  wave  frequency.  Lambda  is  tae 
i 

wavelength,  and  G is  the  gravitational-  acceleration.  The 
wave  surface  is: 


A (x,t) 


= A*cos  (2*Pi<‘x/Lambda+w*t) 
0 i 


= A*  (COS  (w*x/C)  *sin  (w+t)  +sin  (w*x/C)  *cos  (w*t)  ) 
0 i i i i 


(V-7) 

after  substitution  of  Equations  (V-5, 6),  and  the  appropriate 
trigonometric  identities.  The  frequency  experienced  by  the 
craft,  however,  is  the  encounter  frequency: 


L. 
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w =«♦  (C- V*cos  (Beta) ) /C 

e i 


2 

=w-H*V*cos  (Beta)/G 
i i 


(V-8) 

where  V is  the  craft  velocity  and  Beta  is  the  angle  of  craft 
heading  with  respect  to  the  advancing  wave  front,  Beta=0 
degrees  corresponding  to  the  craft  and  wave  front  travelling 
in  the  same  direction  (following  seas)  . 


C.  WAVE  lOfiCING  FUNCTIONS 

Wave  action  introduces  a vertical  force  in  heave  and  a 
moment  about  the  center  of  gravity  in  pitch. 

It  is  assumed  that: 

1.  wave  generated  forces  and  moments  are  due  solely  to  the 
net  additional  hull  volume  submerged  by  the  wave; 

2.  that  the  wave  does  not  exist  within  the  plenum. 

These  assumptions  appear  justified  due  to  the  small 
magnitude  of  waves  considered,  the  large  inertia  of  tne 
craft,  and  the  filtering  action  of  the  bowseal. 

The  waveheight  as  a function  of  the  x-iocation  along  the 
hull  and  time  t may  be  described  by: 

Bvht  (x,t)  =W vht*sin (w *t  + 2*Pi*x/Lambda) 


e 


( V-  9) 


♦cos,(H*t)  ♦sin  (X*w^C) 
e i 


) 


after  substitution  from  equations  (V-5)  and  (V-6)  and  the 
appropriate  trigonometric  identities.  As  sho;#n  by  Fig  15, 
the  differential  volume  of  both  sidewalls  immersed  by  the 
incident  wave  at  time  t is; 

dV(x,t)  = (1/2)  *Wvht  (x,t)  *2*Ms*dx 


and  the  differential  heave  force  at  time  t is; 
dHbuoyw(x,t)  =-Rho^G*dV  (x,  t) 


{ V-10) 


( V-11) 

which  may,  after  appropriate  substitution  from  Equations 
(V-9,10) , be  integrated  analytically  from  stern  to  bow  as; 

Hbuoyw  (t)  =-2*Hho*G*  (C/w)  *Wvht*sin  (w*L/  (2*C) ) *sin  (w^t) 

i i e 

= HBUOrw*sin  (w^t) 
e 


{ V-12) 


The  differential  moment  due  to  the  differential  volume 
submerged  by  the  incident  wave  at  time  t is; 

dPbuoyw  (x,t)  =Hho*G*  (x-Xcg)  ♦dV  (x,t) 


=Rho*G*  (x-Xcg)  ♦Hvht  (x,t)  ♦Ws*dx 


which  integrates  from  stern  to  bow  as; 


( V-13) 


dd 


Pbuoyw  (t) 


= 2*Rho'«‘  Ws*3*W  vht*cos  (w*t)  * 


2 

( (C/wK  *sin  (w*L/(2*C)  ) 
1 i 


- (C*L/(2*w)  ) *cos  (w*L/  (2*C)  ) ) 
i i 


♦ Xcg*Hbuoyw  (t) 


= PBDOYW*co3 (w*t)  +Xcg*HBDOYW*sin (w*t) 
e 9 


( V-14) 


A particular  note  of  caution  should  be  introduced  at  this 

point  with  regard  to  the  integrals  for  wetted  area  and 

wetted  moments  along  the  hull;  the  integrals  of  these 

functions  vanish  for  certain  w , Specifically,  for  neave, 

i 

Equation  (V-12)  vanishes  for  all  time  t when: 


w*L/(2*G)  =n’*'Pi 

i 


1/2 

w = (2*n*?i*S/L)  ,n=1,2,3,... 

i 


( V-15) 


The  pitch  integral  will  likewise  vanish  for  all  time  t if, 
for  a value  of  n such  that  equation  (V-12)  is  zero,  rhe 
transcendental  equation  obtained  by  substitution  of  equation 
(V-15)  into  equation  (V-14),  and  equating  the  latter  to 
zero ; 


2 1/2  2 1/2 
tan  ( (n*Pi*G*L/  (2*C) ) )=  (2*n*Pi*G/  (L*C) ) 


( V-lb) 
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is  also  satisfied.  Thus  the  frequency  response  of  the  craft 
should  exhibit  a certain  ripple  effect  as  these  critical 
frequencies  are  passed. 

The  t«o  previously-defined  wave  integrals  HBUOYW  and 

PBUOYH  make  up  a two-dimensional  B vector  which  couples  the 

wave  force  onto  the  C.G.  acceleration  and  pitch  acceleration 

state  variables,  from  the  two  forcing  functions,  U=sin(w*t) 

1 e 

and  U=cos(w*t).  However,  it  is  convenient  to  redefine  this 
2 e 

relationship  in  terms  of  a single  forcing  function  since  the 

sine  and  the  cosine  are  related  through  the  derivative.  In 
order  to  facilitate  this  relationship,  the  state  matrix  is 
augmented  with  a dummy  state  variable: 

X = d/  (0)  =-w*  (PBUOYW/Iyy)  *U 

10  /dt  2 e 1 


( V-17) 

Thus  the  new  state  x represents,  after  integration,  the 

10 

control  U , multiplied  by  the  appropriate  B (w)  vector  terms 
2 ' e 

of  PBUOYH/Iyj.  Note  that  if  a time  domain  solution  of  wave 

forcing  action  were  desired,  a new  initial  condition  would 
be  required,  representing  the  initial  wave  buoyant  mcment; 
however,  that  does  not  enter  into  the  solution  in  the 
frequency  domain.  The  resulting  state  equations  are  shown 
in  Fig  16.  As  has  been  previously  discussed,  the  frequency 
response  may  now  be  computed  for  any  frequency: 

I -1  I 

G(jw)  =1  1 jw*I-Al  *3(w)  I 

e I e e| 


( V-13) 
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where  G(jw)  is  the  column  vector  of  gains  for  each  state 
e 

variable . 


D.  FBEQaENOi  DOMAIN  SOLUTION 

The  previously-defined  state  equations  were  solved  in 

the  frequency  domain  utilizing  a CSMP  III  computer 

simulation  which  appears  at  the  end  of  the  text.  This 

simulation  was  based  around  the  initialization  portion  of 

the  time  domain  solution  described  in  earlier  chapters; 

however,  rather  than  a call  to  integration  in  the  Dynamic 

portion  cf  the  model,  the  CSilP  III  reserved  variable  TINE 

was  used  only  as  a stepping  variable  to  generate  50 

frequencies  between  w and  w,  the  specified  minimui  and 

0 mar 

maximum  frequencies  respectively,  such  that: 

Delw*TIIlS 

w =w*e 

e 0 

(V-19) 

where  Delw  was  defined  such  that  when  TIME  =1.0,  w=w.  The 

e max 

incident  wave  frequency  then  becomes,  for  ahead  seas,  the 

positive  solution  to  the  quadratic  Equation  (V-8)  , with 

cos  (Beta) =- 1 . The  influence  coefficients  about  the 

steady-state  operating  point  from  the  initialization  portion 

were  copied  onto  a complex  matrix  Ajw,  the  |jwI-A|  matrix, 

e 

and  the  E vector  for  that  particular  incident  wave  frequency 

computed.  The  International  Mathematical  and  Statistical 
Libraries  (IMSL)  subroutine  LEQT2C,  complex  linear 

simultaneous  equatio::  solution,  was  used  to  compute  the 
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vector  of  gains: 

-1 

G(jw)  =1  jwI-A|*B  (w) 

e e e 

(V-20) 

and  the  resultant  absolute  values  of  G(jw)  were  plotted 

e 

against  the  encounter  freguency  w as  decibels  with  respect 

e 

to  the  unit  waveheight  for  heave,  heave  rate,  and  heave 

2 

acceleration,  and  with  respect  to  the  waveslope,  w/G,  for 

i 

pitch,  pitch  rate,  and  pitch  acceleration.  Note  that  chis 

is  a linear  syste:!i;  accordingly,  variations  in  waveheight  do 
not  affect  this  model's  response  characteristics,  as  was 
true  of  the  ncn-linear  Loads  and  Motions  pregram  in  ref  M|. 


E.  COMPAHISCN  WITH  LOADS  AND  MOTIONS  PROGRAM 


Booth  in  ref  |1|  studied  the  frequency  response  of  the  * 
XR-3  Loads  and  Motion  program  in  response  to  seas  developed 
by  the  WAVES  subroutine.  It  is  of  interest  to  compare  the 
linear  system  developed  in  this  study  with  the  results  of 
ref  |1|.  Duplicating  his  work,  three  runs  were  made 
representing  ahead  seas,  cross  seas  and  following  seas  at 
10,  20  and  30  kts  and  a comparison  made  with  his  results. 

Merged  plots  of  pitch  and  heave  were  obtained  for  each  speed 
condition  for  encounter  frequencies  varying  from  1 to  25 
rad/sec.  In  addition,  Cal-Corap  plots  of  pitch  and  heave 
were  obtained  as  functions  of  the  ratio  of  craft  length  to 
wavelength.  Initial  conditions  and  craft  parameters 
selected  for  initialisation  were  identical  to  those  selected 
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for  the  previous  chapter  with  regard  to  time- domain 
validation. 

1 Biiii 

Fig  17  to  Fig  22  show  plots  of  linear  system 
response  in  heave  and  pitch  for  10,  20  and  30  knot  runs  into 
ahead  seas,  alternating  with  plots  reproduced  from  ref  01. 
representing  identical  runs  done  by  Booth  with  the  smallest 
wave  amplitudes.  These  were  selected  as  representing  the 
most  nearly-linear  mode  of  operation  of  the  Loads  and  Motion 
program  in  sea-state.  One  may  observe  the  general 
similarity  of  the  two  simulations,  with  the  speed-dependent 
resonant  peak  of  heave  in  the  linear  model  varying  from  an 
encounter  frequency  of  3.39  rad/sec  at  10  knots  , to  6.45 
rad/sec  at  30  kts,  while  the  non-linear  system  varies  from 
4.0  to  approximately  5.0  rad/sec.  The  frequency  interval 
of  the  non-linear  study  was  considerably  coarser  than  that 
utilized  in  the  linear  model,  1 rad/sec  vs.  less  than  .1 
rad/sec.  The  linear  system~is  significantly  less  responsive 
than  the  Loads  and  Motions  program,  by  11  dB  in  heave  and  6 
dB  in  pitch.  This  is  attributed  to  the  absence  of 
wave-plenum  chamaer  interactions  with  the  plenum  pressure 
Pbbar,  and  indicates  that  rhe  presence  or  absence  of 
significant  wave  action  within  the  plenum  chamber  is  an  area 
of  considerable  impact  on  the  response  of  the  craft,  and 
should  be  included  in  future  studies. 

The  linear  system  is  quite  sharper  in  its  resonant 
peak  in  pitch,  as  should  be  expected  from  the  significantly 
reduced  pitch  damping  previously  demonstrated  by  the  linear 
model  in  the  time  domain  validation.  This  accounts  in  part 
for  the  tendency  of  the  linear  model  to  decrease  in  response 
slightly  from  an  encounter  frequency  of  1 rad/sec,  until 
just  before  the  onset  of  resonance,  unlike  the  Loads  and 


93 


[ 

1 Motions  data;  the  half-power  points  in  the  Loads  and  Motions 

I frequency  resonance  are  more  widely  separated  than  are  those 

of  the  linear  system,  and  this  in  part  also  accounts  for 
this  difference. 

The  linear  system  shows  the  distinctly  periodic 
frequency  response  predicted,  with  a sharp  null  occurring  at 

incident  wave  frequencies,  w,  of  3.14,  4.4,  and  5.38 

1 

rad/sec.  These  frequencies  correspond  to  a value  of  n,  in 

eqn  (V-15),  cf  1,  2 and  3.  These  nulls  are  quite  sharp,  and 
are  not  as  well  defined  in  the  Loads  and  Motions  studies  due 
to  the  coarser  frequency  resolution,  but  are  nonetheless 
visible  in  all  runs. 

The  Loads  and  Motions  studies  also  demonstrate  an 
interesting  form  of  non-linear  response  with  regard  to 
resonant  pitch  frequency.  For  small  wave  amplitudes,  less 
than  .1  ft,  the  Loads  and  Motions  resonant  pitch  frequency 
varies  with  craft  speed  in  nearly  an  identical  manner  with 
the  linear  system,  with  the  slight  exception  of  the  30  kt 
study;  however,  the  only  available  data  for  that  run  was  at 
a wave  amplitude  of  .2  ft.  For  wave  amplitudes  greater  than 
.1  ft,  the  resonant  peak  occurs  at  a frequency  of  4 rad/sec, 
regardless  of  craft  speed.  A tentative  explanation  for  this 
behavior  is  herein  advanced:  as  the  craft  resonant 

frequency  is  approached,  bowseal  gapping  commences  for  the 
larger  wave  amplitudes,  with  attendant  loss  of  plenum 
airmass;  the  loss  of  seal  force  and  plenum  pressure  cause 
the  craft  to  settle  back  into  the  water,  but  a fixed  amount 
of  time  must  elapse  before  the  fans  can  replenish  the  plenum 
airmass  and  restore  normal  system  operation.  Thus  the  4 
rad/sec  period  observed  for  larger  wave  amplitudes 
represents  the  fastest  time  (.66  sec)  over  which  gapping 
can  occur  and  the  fans  sufficiently  replenish  the  airmass  to 
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perait  it  to  recur.  This  non-linear  response  is  soaewhat 
analogous  tc  a relaxation  oscillator,  and  represents  an 
interesting  area  for  future  study,  as  such  a non-linearity 
might  easily  be  modeled  in  the  Kyguist  plane,  utilizing  the 
modified  Nycuist  stability  criteria  to  predict  this 
non-linear  behavior. 

a.  10  Kts  Ahead  Seas 

Figs  17  and  18  show  the  superimposed  heave  and 
pitch  response  of  the  linear  and  Loads  and  Motions  program, 
respectively,  for  an  ahead-seas  run  at  a craft  velocity  of 
10  kts.  The  linear  model  represents  the  response  of  the 
craft  in  heave  and  pitch  at  50  discrete  frequencies  ranging 
exponentially  from  1.0  to  25.0  rad/sec.  The  non-linear  data 
reproduced  from  ref  j^lj  represents  8 discrete  frequencies 
over  that  same  range  at  1.0,  2.0,  3.0,  4.0,  5.0,  6.0,  8.0, 
and  10.0  rad/sec,  at  a waveheight  of  .2  ft.  Although  this 
waveheight  was  the  smallest  used  by  Booth  for  runs  at  these 
speeds,  it  is  still  sufficiently  large  to  excite  the 
aforementioned  non-linear  response.  This  behavior  makes 
itself  manifest  as  the  4.0  rad/sec  resonance  in  pitch, 
compared  iith  3.39  rad/sec  for  the  linear  model.  The  Loads 
and  Motions  program  shows,  for  small  weight  disturbances  at 
these  speeds,  a natural  pitch  period  of  1.64  secs,  or  3.7 
rad/sec.  As  previously  noted,  the  magnitude  of  the  linear 
system  pitch  response  agrees  fairly  well,  within  3 dE,  of 
the  non-linear  Loads  and  Motions  data,  with  a sharper  peak 
due  to  the  lighter  linear  damping;  the  heave  magnitudes 
differ  considerably,  by  -15  dB,  due  to  the  absence  of  plenum 
waves  in  the  linear  system. 

liith  the  exception  of  the  heave  magnitude 
difference  and  the  shift  in  resonance  due  to  the  non-linear 
behavior,  tie  shape  of  the  heave  response  curves  of  the  two 


95 


AD-A044  325  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CALIF  F/6  1/3 

investigation  of  a linear*  TWO-DEGREE-OF-FrEEDOM  SIMULATION  OF  — ETCCU) 
JUN  77  L F MCINTYRE 

UNCLASSIFIED  NL 


ACVV0443^ 

a 

r 



i. 

> 

1* 

- 

a 

:: 

a 

r 

PI 

— i 

□ 

□ 

gp 

a 

a 

a 

1 

[ 

1 

: 

\ 1 

-- — 

< 

. i 

1 

1 

' f 

1 

i 

_..i 

f 1 

t 

— »-»T- •* 

.■ . T 

f 



J 

’1 

t*- 

■‘i  i i .i- 

i ' • . . ’ 

' 1 

■fe:::! 

1 

L_  . 

END 

DATE 

fILHCO 

10-77 

DDC 

I 


models  ccmpare  quite  favorably.  A shallow  null,  of 
undetermined  origin:,  appears  in  both  heave  curves  at  3.0 
rad/sec,  just  slightly  before  the  onset  of  resonance.  A 
very  sharp  null  occurs  at  8.36  rad/sec  in  both  models;  this 
is  at  an  incident  wa»e  frequency  of  3.15  rad/sec,  or  a 
wavelength  of  20  ft.  As  predicted  by  Equations  (V-15,16), 
further  nulls  occur  at  a wave  incident  frequency  of  4.46  and 
5.58  rad/sec,  corresponding  to  wavelengths  of  10  and  6.6 
feet  respectively,  but  the  Loads  and  Motion  encounter 
frequencies  did  not  range  sufficiently  high  to  observe  them. 


Pitch  also  shows  general  agreement;  but  Booth’s 
Loads  and  Motion  data  shows  a null  at  6.0  rad/sec,  vs.  11.54 
rad/sec  for  the  linear  system.  This  null  may,  however,  be 
related  to  the  non-linear  response  previously  described, 
representing  a frequency  at  which  the  phase  relationship 
between  plenum  airmass  replenishment  and  pitch  buoyancy  due 
to  waves  is  such  that  cancellation  occurs.  It  is 
suspiciously  located  at  a 3/2  multiple  of  the  non-linear 
response  frequency  of  4.0  rad/sec.  In  addition.  Booth  noted 
this  identical  discrepancy  in  comparison  with  the  linearized 
Aerojet  XE-3  model;  the  Aerojet  model  also  predicted  a null 
in  the  vicinity  of  12  rad/sec,  which  Booth  was  unable  to 
observe  with  the  Loads  and  Motions  program,  but  which  is 
distinctly  visible  in  the  linear  model  at  11.54  rad/sec. 
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b.  20  Kts  Ahead  Seas 


Figs  19  and  20  show  the  linear  and  Loads  and 
Motions  heave  and  pitch  frequency  response  for  a 20  let  run 
into  ahead  seas.  The  Loads  and  Motions  data  is  for  a wave 
height  of  .1  ft  and  is  more  nearly  linear  in  behavior  than 
previously  observed  at  10  kts.  The  range  of  frequencies 
spanned  by  the  linear  system  is  the  same  as  before;  the 
Loads  and  Motion  data  includes  two  additional  points  at  15 


r 


and  20  rad/sec  in  heave  response.  The  linear  pitch  response 
is  snaller  than  that  observed  at  10  kts,  being  about  6 dB 
les.<;  than  the  non-linear  data,  vhile  the  heave  response 
remains  consistently  less  than  the  non-linear  data  by  12  dB 
as  at  10  kts. 


The  shape  of  the  two  sets  of  frequency  response 
data  is  in  quite  good  agreement  for  these  conditions,  as  the 
Loads  and  Motions  is  apparently  in  a nearly-linear  mode  of 
operation.  The  pitch  data  shows  a distinct  flattening  at 
the  peak  of  resonance  in  both  models,  although  this  occurs 
somewhat  earlier  in  the  non-linear  data  (4.0  rad/sec)  than 
in  the  linear  data(4.68  rad/sec),  due  to  the  more 
sharply-peaked  pitch  response  of  the  latter.  Both  models 
resonate  at  5.0  rad/sec,  and  fall  off  above  resonance  with 
identical  roll-offs  of  12  dB  per  octave,  after  the  initially 
steeper  roll-off  of  the  linear  system  from  peak  resonance. 
Heave  shows  the  characteristic  null  in  both  models  at  an 
encounter  frequency  of  14  rad/sec;  this  is  a wave  frequency 
of  3.  15  rad/sec,  or  a wavelength  of  20  ft.  The  linear 
system  predicts  a pitch  null  at  20.0  rad/sec,  but  the  Loads 
and  Motions  data  does  not  extend  sufficiently  high  to 
observe  this. 


c.  30  Kts  Ahead  Seas 


figs  21  and  22  show  the  linear  and  non-linear 
frequency  response  for  a 30  kt  ahead-seas  run.  As  in  the  10 
kt  run,  the  smallest  wave  amplitude  utilized  in  Booth’s 
study  was  .2  ft;  however,  it  is  considered  that  the 
nose-down  planing  moment  which  is  quite  strong  at  high 
speeds  restrains  bowseal  gapping,  and  the  Loads  and  Motions 
data  represents  conditions  for  the  non-linear  response  which 
is  just  on  the  verge  of  becoming  established.  This  is 
supported  by  the  observation  that  pitch  resonance  occurs 
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between  the  frequency  previously  observed  as  the 
characteristic  non-linear  response,  and  linear  resonance, 
between  4.0  and  5'.  0 rad/sec;  furthermore,  the  heave 
resonance  is  weaJcly  peaiced  at  a slightly  higher  frequency 
than  pitch  response,  approximately  5.0  rad/sec,  still  closer 
to  the  predicted  linear  resonance  at  6.46  rad/sec. 

The  linear  model  shows  still  lower  magnitudes  of 
pitch  response,  6 to  8 dB  below  the  Loads  and  Motion  data. 
The  heave  magnitudes  remain  consistently  12  dB  below  the 
Loads  and  Motion  data. 

Pitch  resonates  at  6.46  rad/sec  in  the  linear 
system,  and  between  4.0  and  5.0  rad/sec  in  the  non-linear 
Loads  and  Motions  data.  Roll-off  in  pitch  remains  the  same 
above  resonance,  approximately  12  dB  per  octave  in  both 
models,  but  heave  response  in  both  models  falls  off  more 
shallowly,  approximately  10  dB  per  octave  in  both  models. 
The  heave  null  due  tp  the  20  ft  wavelength  is  still  apparent 
in  the  linear  model,  representing  an  encounter  frequency  of 
20.0  rad/sec,  but  the  loads  and  Motions  data  did  not  extend 
sufficiently  high  in  frequency  for  this  null  to  be  observed. 


d.  Frequency  Response  vs.  Wavelength 


Pigs  23  and  24  are  plots  of  heave  response,  vs. 
the  ratio  cf  craft  length  to  wavelength  for  the  linear 
system,  and  vs.  the  ratio  of  wavelength  to  craft  length  for 
the  Loads  and  Motion  data  reproduced  from  ref  Both 
plots  are  superimposed  with  the  data  from  all  three  velocity 
runs  of  10,  20,  and  30  Jets.  The  difference  in  scale 
selected  for  plotting  in  the  linearized  model  precludes 
crowding  at  the  low  end  of  the  x-axis,  since  the  bulk  cf  the 
frequencies  studied,  3.15  rad/sec  and  above,  are  shorter 
than  craft  length. 


98 


A 


The  linear  system  heave  response  shows  distinct 
nulls  for  waves  whicii  are  integer  multiples  or  submultiples 
of  the  craft  length,  and  may  be  observed  at  L/Lambda= 
1/4, 1.0, 2.0,  and  3.0;  less  distinct  nodes  may  be  observed  at 
L/Lambda=  3/2,  5/2.,  and  7/2.  The  corresponding  nulls  at 
Lambda/L=4 . 0, 1 . 0,  and  .5  are  visible  in  Fig  24,  but  scale 
crowding  and  the  coarser  freguency  resolution  utilized  by 
Booth  render  these  less  obvious.  The  sharp  node  at 
L/Lambda=.3  corresponding  to  a wave  frequency  of  1.75 
rad/seC  is  net  the  result  of  the  periodic  hull  integrals; 
this  wave  freguency  corresponds  to  the  encounter  freguency 
of  3.4  to  6.C  rad/sec,  and  is  the  craft's  natural  pitch 
freguency.  It  may  be  observed,  also  as  a node,  in  the  pitch 
response.  The  natural  pitch  freguency  variation  with  speed 
is  obscured  in  this  plot,  since  as  the  speed  increases,  the 
natural  pitch  freguency  increases;  but  this  is  an  encounter 
frequency,  and  the  wave  incident  freguency  which  produces 
this  encounter  freguency  decreases  with  increasing  speed. 
Therefore,  the  two  effects  nearly  cancel. 

Figs  25  and  26  show  pitch  response  for  the 
linear  system  and  the  non-linear  loads  and  Motion  data  of 
ref  W.  plotted  as  was  heave.  Note  that  the  nulls  and 
nodes  are  roughly  18Q0  out  of  phase,  a node  in  pitch 
corresponding  to  null  in  heave  and  vice-versa. 


2.  Cross-seas  Runs 


Figs  27  to  31  are  presentations  of  heave  and  pitch 
freguency  response  of  the  linear  system  and  the  Loads  and 


Notions  data  from  ref 


f'] 


at  20  and  30  kts  under  cross-seas 


conditions 


(w=y).  • 

e 1 


Comments  made  regarding  the  the 
ahead-seas  runs  are  generally  applicable  here,  but  the 
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comparison  of  the  two  systems  is  in  no  manner  nearly  as 
striking  as  thar  for  the  ahead-seas  runs.  This  is 
attributed  to  the  significant  effects  of  roll  and  yaw  modes 
on  the  craft  dynamics  under  these  conditions,  which  are 
completely  neglected  in  this  linear  model. 

3 • Following  Seas  Runs 


•j 

I 

i( 

li 
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Booth  acknowledged  little  agreement  between  his  work 
with  the  Loads  and  ilotion  program  and  the  Aerojet  XR-3 
linear  model  frequency  response  curves;  this  is  apparently 
due  to  an  error  in  calculating  the  encounter  frequency. 
Higher  frequency  waves  propagate  at  a slower  velocity,  so 
that  for  a craft  traveling  at  a given  velocity,  as  the 

incident  wave  frequency  w increases,  a frequency  will  be 

i 

attained  above  which  the  craft  will  actually  outrun  the 

following  wave,  and  plow  into  such  a wave  bows-on.  Thus  the 
craft  will  actually  experience  an  ahead-sea  condition  above 

that  wave  frequency.  The  w corresponding  to  the  maximum 

i 

positive  w that  the  craft  may  experience  for  following  seas 
(cos  (Beta)  = 1 . 0)  is  given  by: 

2 

d/dw(  w-V*w/G)  =0.0 
1 i i 


(V-21) 


I 
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= G/  (2*V) 


w 


1 


Substituting  Equation  (V-22)  into  (V-8)  yields: 

2 

w =G/ (2’»V)- (V/G)  ♦ (G/ (2»V)  ) 

emax 

=G/  (U*V) 


(7-22) 


(V-23) 

Incident  wave  frequencies  above  or  below  the  w corresponding 

i 

to  w result  in  a lowar  w,  and  if  w is  increased 

emax  e i 

suf f icie.^tly , past  w=G/V,  w will  become  negative,  indicatinq 

i e 

that  the  craft  is  overrunning  the  wave.  Accordingly,  a great 
deal  of  judiciousness  must  be  utilized  when  interpreting  the 
results  of  a following  seas  run,  when  plotted  against 
encounter  frequency  w. 


Accordingly,  the  following-seas  frequency  spectrum 

was  generated  in  a manner  different  from  that  utilized  on 

previous  runs.  M was  varied  exponentially  over  the  range  of 

i 

1 to  25  rad/sec,  and  the  corresponding  w was  calculated. 

e 

Hence,  the  results  are  not  comparable  directly  with  the 

results  of  Booth,  but  are  presented  in  Fig  32  to  Fig  33  for 
information  only,  at  v=20  lets. 
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Figure  16  - STAJE  VARIABLE  REPRESENTAIIOy  OF  THE  WAVE 

FORCING  FUNCTION 
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LINEAR  SYSTEM  FHEgUENCY  RESPONSE, 
10  KTS,  AHEAD  SEAS 
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I 

i 

V=10  kts,  A=0.2  ft,  i 

Ahead  Seas 

0 = Heave  Frequency  Response 
O = Pitch  Frequency  Response 

Figure  18  - LOADS  AND  MOTION  FE2QUENCY  RESPONSE, 

10  KTS,  AHEAD  SEAS 
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LINEAR  SYSTEM  FREQUENCY  RESPONSE, 
20  KIS,  AHEAD  SEAS 
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V=20  kts,  A=0.1  ft, 

Ahead  Seas 

0 = Heave  Frequency  Response 
0 = Pitch  Frequency  Response 

Figure  20  - LOADS  AND  MOTION  FEE^UENCY  RESPONSE, 
20  KTS,  AHEAD  SEAS 
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figure  21  - LINEAB  SYSTEM  FREQUENCY  RESPONSE, 
30  KTS,  AHEAD  SEAS 
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V=30  kts,  A=0.2  ft, 

Ahead  Seas 

0 = Heave  Frequency  Response 
O = Pitch  Frequency  Response 

Figure  22  - LOADS  AMD  MOTION  FREQUENCY  RESPONSE, 
30  KTS,  AHEAD  SEAS 
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Figure  23  - LINEAH  HEAVE  FREQUENCY  RESPOHSE  VS.  L/LAHDDA 
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Figure  27 


LINEAR  5YSTE.1  FSEiJUENCY  RESPONSE, 
20  Krs,  CROSS  SEAS 
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Figure  28  - LOADS  AND  MOTIONS  HEAVE  RESPONSE, 
20  AND  30  KTS,  CROSS  SEAS 
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V=20  lets,  A=0.2  ft, 

Abeaun  Seas 

O = Pitch  Frequency  Response 
A = Roll  Frequency  Response 

Piyurt  29  - LOADS  AND  MOTION  PITCH  RESPONSE, 
20  KTS,  CROSS  SEAS 
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Figure  30  - LINEAR  SYSTEM  FREQUENCY  RESPONSE, 
30  KTS,  CROSS  SEAS 
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LINEAR  SYSTEM  FREQUENCY  RESPONSE 
20  KTS,  FOLLOWING  SEAS 
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V=20  kts,  A=0.2  ft 
Following  Seas 

0 = Heave  Frequency  Response 
O = Pitch  Frequency  Response 

Figure  33  - LOADS  AND  MOTION  FREQUENCY  RESPONSE, 
20  KTS,  FOLLOWING  SEAS 
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VI 


PREaUENCY  DOMAIN  RSDDCTION  OP  XE^l  DATA 


This  chapter  develops  the  techniques  for  fast  efficient 
reduction  of  XH-3  testcraft  data  froQ  the  time  domain  to  the 
frequency  domain.  Some  more  prominent  noise  sources  are 
identified,  and  the  mode  of  noise  transmission  is 
hypothesized.  A very  limited  amount  of  pitch  data  is 
collected  and  closely  examined  using  these  techniques  and 
compared  with  the  theoretical  development  of  the  preceding 
chapters . 


A.  DATA  ACQCISITION  SYSTEM 


The  data  acquisition  system  consists  of  the  standard 
semi“per lanently  installed  Pemco  tape  recorder.  This  is  a 
14-channel  PM  tape  system,  with  output  0-1  volt  positive.  A 
19-channel  patch  panel  in  the  tape  recorder  compartment, 
located  just  aft  of  the  cockpit,  permits  selection  or  tne 
data  to  be  recorded; voice  edge  track  recording  directly  from 
the  craft's  radio/ICS  system  permits  relatively  simple 
documentation  of  experimental  runs  to  be  maintained. 


B.  FREQUENCY  DOMAIN  DATA  HEDUCIICN  SYSTEM 


The  Eio- Engineering  Labo 
Department,  NPS,  Monterey, 
for  spectral  analysis  of 
consisting  cf  a PDP-11 


ratory.  Electrical  Engineering 
CA,  has  a very  powerful  system 
analog  data.  This  system. 
Mini-computer  and  Time-Data 
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Corporation  Superpanel-50  Sp.ectrum  Analyzer,  was  utilized 
for  time-to-f requency  domain  data  reduction  of  the  analog 
tape  data.  Ibis  system  is  capable  of  performing  a fast 
Fourier  Transform  (FFT)  on  analog  data  of  magnitudes  of  +5'/, 
from  a maiimum  upper  frequency  of  10  Hz  to  10  MHz,  for  up  to 
80S4  discrete  time  samples.  However,  caution  must  be 
exercised  to  avoid  memory  overflow  and/or  program  interrupt 
of  the  EDP-11  for  samples  in  excess  cf  2048.  In  addition, 
this  entire  process  may  t>e  repeated  up  to  2048  times,  and 
the  resultant  Fourier  coefficients  averaged,  for  further 
reduction  cf  transient  frequency  components.  The  resultant 
transformed  data  may  be  output  as  Fourier  coefficients, 
power  spectral  coefficients,  and  as  an  amplitude  histogram; 
output  devices  include  a storage  C3T  and  a Hewlett-Packard 
70C4-B  X-y  plotter  which  will  reproduce  the  display  on  the 
CRT  on  ccmmand  from  the  panel.  In  addition  the  time  domain 
data  being  analyzed  may  be  observed  on  a time-base  storage 
oscilloscope,  and  monitored  throughout  the  transformation 
process.  Processing  times  vary,  of  course,  with  upper 
frequency,  cumber  of  samples,  and  number  of  averaging 
"windows",  fiom  less  than  a second  to  several  minutes.  For 
further  information  see  ref  ^8^ , the  operating  instructions 
for  the  Superpanel-50.  Fig  34  shows  the  details  of  the 
setup  for  data  reduction  to  the  frequency  domain. 


C.  FREQUEUCY  DOMAIN  ANALYSIS 


XH-3  data  was  collected  on  an  experimental  run  conducted 
20  January,  1977  at  Lake  San  Antonio  for  the  purpose  of 
performing  frequency  domain  analysis,  using  the 
previously-described  FFT  techniques.  Two  analyses  were 
performed ; 

1.  a broad-band  noise  analysis,  examining  the  frequency 
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spectra  of  all  channels  up  to  200  Hz,  as  various 
suspected  noise  sources  were  brought  up  and  shut  down, 
or  varied  in  power  setting;  the  purpose  of  this 
experiment  was  to  identify  offending  noise  sources  by 
their  characteristic  signature  in  the  frequency  domain; 

2.  a Icw-f requency  analysis  of  craft  characteristics  in 
pitch  and  pitch  rate,  for  15  let  runs  into,  across,  and 
with  the  prevailing  wave  directions.  The  frequency 
spectra  of  these  channels  were  analyzed  from  .025  to 
25  Hz,  for  the  purpose  cf  comparison  with  ref  I'J  and 
the  werk  cf  the  preceding  chapters. 

^ Analysis 

Heise  data  from  the  XR-3  was  analyzed  over  a 
frequency  range  up  to  200  Hz  and  512  samples,  averaged  over 
8 windows;  this  gave  a total  sample  time  of  approximately  30 
seconds.  A Butterworth  filter  was  used  to  reject  signals 
above  200  Hz  to  prevent  aliasing  of  higher  frequency 
components  into  the  lower  range  of  interest  due  to 
discretization.  In  addition,  a battery/ri^eostat  combination 
was  used  to  bias  out  the  1/2  volt  offset  of  the  F3  tape, 
since  , as  previously  noted,  the  FM  tape  output  range  is  0-1 
volt  positive,  while  the  input  to  the  Analog- to-Digital 
converter  of  the  Superpanel-50  is  ^5V;  this  consideranly 
improves  amplitude  resolution  of  the  transformed  data. 

A representative  output  of  frequency  spectra  appears 
in  Fig  35;  the  remainder  may  be  examined  in  Appendix  A.  The 
data  from  all  channels  are  basically  identical  in  nature, 
and  too  repetitive  to  include  in  the  main  text.  As  lay  be 
expected,  the  craft  main  engines  and  the  60-Hz  power  supply 
contribute  the  major  share  of  noise.  As  power  is  varied,  the 
main  engine  noise  may  be  observed  to  vary  from  20  Hz  tc  100 
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Hz,  while  the  power  supply  contributes  significant  noise  a-*- 
55-60  Hz  and  its  harmonics.  Following  these  sources  in 
order  of  magnitude  are  the  lift  fans,  which  may  be  observed 
in  the  vicinity  of  40  Hz. 

An  attempt  was  made  to  ascertain  the  mode  of 
transmission  of  this  noise,  specifically,  whether  by 
electrical  means,  or  by  mechanical  vibration  of  the  sensor 
package.  Accordingly,  the  gyro  package  containing  pitch, 
roll,  and  yaw,  and  pitch  rate,  roll  rate,  and  yaw  rate  gyros 
was  removed  and  placed  on  a Ling  shaker  table,  and  exposed 
to  vibrations  of  1/2  to  1 G,  over  a frequency  range  of  5-200 
Hz.  The  RMS  shaker  table  displacement  was  measured  with  a 
calibrated  Bentley  accelerometer,  and  compared  with  the  RMS 
voltage  output  of  the  individual  gyros.  It  was  hoped  to 
obtain  frequency  response  curves  of  the  gyro  packages,  but 
they  seemed  quite  insensitive  to  vibrations  in  this 
frequency  range,  and  the  frequency  curves  showed  significant 
scatter  due  to  the  low  signa 1-to- noise  ratio.  They  were  not 
considered  worthy  of  inclusion. 

It  was  concluded,  on  the  basis  of  this  experiment, 
that  mechanical  vibration  could  not  transmit  such 
overwhelming  amounts  of  noise  as  was  observed  being  injected 
into  the  XR-3  data  collection  system.  The  medium  of 
transmission  is  thus  electrical;  this  conclusion  is  further 
supported  by  the  observation  that  the  eight  engines 
installed  aboard  the  XR-3  have  no  ignition  noise 
suppressors;  likewise,  the  60  Hz  power  supply  is  unfilterei. 
Furthermore,  while  bonding  straps  run  the  length  of  the 
craft,  the  possibility  of  extensive  AC  ground  loops  is 
always  present  in  a non-conducting  hull,  such  as  the  XR-3. 
These  areas  should  receive  some  close  scrutiny  in  the 
future,  for  a tenfold  reduction  in  noise  injection  into  th<= 
data  collection  system  should  be  possible,  by  means  of 
ignition  noise  suppression,  and  the  installation  of  .1  uF 
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capacitors  at  each  AC  outlet,  at  a relatively  insignificant 
cost  in  time  and  money. 

2 • lov-trcquencY  Craft  Dyn^gics 

Ihree  runs  were  made  during  the  aforementioned  data 
collection  rrns,  in  which  pitch  and  pitch  rate  were  recorded 
for  15  kt  runs  into,  across,  and  with  the  prevailing  seas. 
Qualitative  sea  conditions  at  the  time  of  observation  were 
noted  as  wave  heights  of  3 to  6 inches,  with  a wavelength  of 
approximately  6 feet;  this  wavelength  represents  an  incident 
wave  frequency  of  approximately  2.2  rad/sec,  corresponding 
to  encounter  freguencies  of  6.2,  2.2,  and  -1.8  rad/sec,  for 
ahead,  cress,  and  following  seas  at  15  kts.  Data  was 
reduced  in  a manner  similar  to  that  used  for  wideband  noise, 
with  the  exception  tbat  a maximum  frequency  of  25  Hz,  1024 
samples  with  no  averaging  was  utilized.  > A difficulty 
developed  in  that  a Butterworth  filter  with  cut-off  below 
100  Uz  was  not  readily  available  at  the  time.  However,  when 
cne  is  constrained  from  lowering  the  bridge,  it  may 
nevertheless  be  possible  to  raise  the  river.  In  this 
particular  case,  the  tape  was  run  at  four  times  recorded 
speed,  or  1 1/2  inches  per  second.  The  FFT  was  set  for  100 
Hz,  as  was  the  Butterworth  filter,  and  the  desired 
resolution  of  1024  samples  selected.  This  procedure  had  the 
beneficial  side-effect  of  reducing  the  real  elapsed 
processing  tine  of  the  Super panel-50/PDP- 1 1 combination  by  a 
factor  of  four.  The  resultant  transformed  data  was 
scale-expanded  after  computation,  so  that  only  the  portion 
frem  0-5  Hz  (0-1.25  actual  Hz)  was  plotted. 

These  admittedly  unorthodox  and  awkward  procedures 
were  necessitated,  not  only  by  the  lack  of  suitable 
filtering,  but  also  by  difficulties  experienced  in  obtaining 
proper  operation  of  the  PDP-11  for  very  low  frequencies  at 
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that  particular  time.  Long  sample  times  combined  with 
frequent  inexplicable  program  interrupts  rendered  this 
particular  session  most  frustrating.  However,  ventilation 
and  overheating  problems  were  apparently  the  major  source  of 
difficulties  experienced  with  this  particular  exercise,  and 
future  data  reductions  in  the  frequency  domain  should 
require  none  of  these  awkward  circumlocutions. 

a.  Pitch  and  Pitch  Hate-Ahead  Seas  Hun 


figs  36  and  37,  the  reduced  FFT  of  pitch  and 
pitch  rate  for  the  15  kt  ahead-seas  run,  may  be  compared 
with  Figs  3ti  and  39,  the  linear  system  theoretical  response 
in  pitch  and  pitch  rate.  The  linear  system  was  plotted  as 
dB  with  respect  to  waveslope  logarithmically  against 
encounter  frequency,  on  a range  of  .1  to  7.85  rad/sec.  The 
observ-ed  data  are  relative  magnitudes  plotted  linearly  vs. 
encounter  frequency,  on  a range  of  0 to  7.85  rad/sec.  The 
unusual  scale  was  necessitated  for  the  observed  data,  to 
obtain  the  conversion  to  rad/sec;  in  Hz,  the  scale  is  .125 
Hz/in.  A qualitative  comparison  of  the  two  plots  shows  a 
rather  remarkable  similarity;  the  FFT  of  pitch  shews  a 
falling  off  in  response  from  the  very  low  frequency  end  to  a 
local  minima  at  2.0  rad/sec  with  the  onset  of  resonance 
occurring  in  the  vicinity  of  2.75  rad/sec;  resonance  peaks 
at  3.9  rad/sec  and  is  completed  at  approximately  5.0 
rad/sec.  Two  small  nodes  may  be  observed  at  2.36  and  6.28 
rad/sec,  the  latter  corresponding  to  the  encounter  frequency 
of  the  qualitatively-observed  wavelength  at  the  time  of  this 
run;  this  was,  as  previously  noted,  an  incident  wave 
frequency  of  2.2  rad/sec. 


The  linear  model  replicates  the  coserved  data  in 
pitch  almost  precisely,  with  the  exception  of  the  two  nodes 
at  2.2  and  b.28  rad/sec,  due  to  the  dominant  local  wave 
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amplitudes  in  the  observed  data. 


As  may  be  expected,  pitch  rate  data  is  somewhat 
more  noisy  than  pitch  data,  although  not  excessively  so. 
observed  pitch  rate  shows  an  unusual  double  peak  in 
resonance,  which  also  occurs  in  the  other  two 

sea-conditions.  The  first  peak  occurs  at  h.35  rad/sec, 

slightly  higher  in  frequency  than  the  observed  pitch  data, 
while  the  second  peak  occurs  at  5.25  rad/sec  and  is  of 
unexplained  origin.  This  peak  may  be  indicative  of  the 
onset  of  the  non-linear  behavior  discussed  in  the  preceding 
chapter,  but  that  is  only  a speculative  observation. 

The  linear  system  pitch  rate  data  shews  a 
sharper  resonance  than  did  linear  pitch,  the  former  being 
markedly  peaked  at  4.26  rad/sec,  while  the  latter  was 

relatively  flat  from  3.9  to  4.26  rad/sec.  The  agreement  in 
the  sharpness  of  pitch  rate  response  between  the  observed 
data  and  the  linear  system  is  quite  good.  Resonant  rise 
also  occurs  earlier  in  the  linear  model,  at  approximately 
2.0  rad/sec;  but  there  is  no  indication  of  the  peak  at  5.25 
rad/sec  which  occurs  in  the  observed  data. 

Ho  significant  nodes  or  nulls  occur  xn  this 

frequency  range;  the  pitch  node  at  n=1(  w=3.15  rad/sec)  does 

i 

not  occur  until  an  encounter  frequency  of  9.75  rad/sec. 

Overall,  the  observed  and  linear  data  ccapare 
remarkably  well  under  these  conditions. 


b.  Pitch  and  Pitch  Rate-Cross  Seas  Run 

Despite  the  fact  that  cross-seas  runs  induce 
modes  of  motion  in  roll  and  yaw  that  are  unaccounted  for  in 
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this  linear  model,  it  was  considered  of  interest  to  compare 

observed  cross-seas  data.  Figs  40  and  41,  with  Figs  42  and 

43,  the  linear  data  for  w=w. 

e i 

The  observed  pitch  FFT,  Fig  40,  is  more 
distinctly  resonant  at  4.0  rad/sec  under  these  saa 

conditions  than  in  the  ahead-sea  runjresponse  decreases 
continuously  from  the  very  low  frequency  end  until  almost 
the  onset  cf  resonance  at  approximately  3.5  rad/sec. 
Resonance  occurs  at  4.0  rad/sec,  and  a distinct  null  occurs 
at  5.25  rad/sec.  The  observed  predominant  wavelength,  which 
as  noted  previously  corresponds  to  an  incident  wave 
frequency  cf  2.2  rad/sec,  is  distinctly  visible  at  that 
frequency . 

Comparison  with  the  linear  system  pitch  response 
curves  cf  Fig  42  is  as  striking  as  the  ahead-seas  runs,  the 
more  so  because  of  the  modes  of  motion  induced  in  the  craft 
which  are  not  in  the  linear  model.  As  in  the  FFT  data, 
linear  pitch  response  decreases  until  almost  the  onset  of 
resonance,  with  a pronounced  null  at  3.5  rad/sec.  This  null 
does  not  appear  in  the  observed  pitch  data,  but  is 
distinctly  visible  in  the  observed  pitch  rate  data.  linear 
system  resonance  is  sharper  under  these  sea  conditions,  as 
is  the  observed  data,  and  occurs  at  4.26  rad/sec.  Two 
high-frequency  nulls  at  5.25  and  7.0  rad/sec,  and  a node  at 
6.0  rad/sec  , are  distinctly  visible  in  the  observed  data, 
although  less  prominently. 

The  observed  pitch  rate  data  shows  a sharper 
rcll-off  from  the  low  frequency  end  of  the  spectrum,  and  a 
flatter  response  thereafter  until  the  onset  of  resonance, 
than  does  pitch.  The  distinct  double  peak  previously  noted 
in  the  ahead-sea  run,  is  distinctly  visible  in  this  run 
also. 


128 


r 


Ccmparispn  of  the  pitch  rate  FFT,  Fig  41,  with 
the  linear  system  pitch  rate  response  of  Fig  43  shows  the 
same  good  agreement  as  did  the  pitch  data;  the  low  frequency 
roll-off  characteristics  are  very  much  the  same,  and  the 
null  predicted  at  3.5  rad/sec  by  the  linear  model  is 
prominently  visible  in  the  FFT.  However,  the  higher 
frequency  nodes  and  nulls  are  much  less  distinct. 

c.  Pitch  and  Pitch  Hate-Following  Seas 

Due  to  the  difficulties  previously  experienced 
in  interpreting  results  of  following  seas  runs  with  3coth*s 
woric,  only  a most  general  comparison  of  observed  pitch  and 
pitch  rate  FFT's,  shown  in  Figs  44  and  45,  and  the  linear 
system  responses,  shewn  in  Figs  46  and  47,  will  be  made. 
Hote  that  due  to  the  technique  utilized  in  generating 
encounter  frequencies  in  the  linear  model,  the  range  of 
frequencies  plotted  in  the  linear  model  is  much  greater,  (.3 
rad/sec  to  -40  rad/sec) , resulting  in  a much  coarser 
frequency  resolution  in  the  range  of  interest.  Negative 
encounter  frequencies,  as  noted  in  the  preceding  chapter, 
represent  a modified  ahead-sea  condition. 

The  linear  system  predicts  a sharp  roil-eff  in 
pitch  response  from  the  very  low  frequency  end  to  about  -.25 
rad/sec,  and  a relatively  flat  response  from  there  until  tae 
onset  of  resonance,  at  about  -2.0  rad/sec.  Resonance  occurs 
quite  sharply  at  -4.11  rad/sec.  No  nodes  or  nulls  in  this 
frequency  range  are  predicted  by  the  linear  model.  The 
observed  pitch  data  is  in  general  agreement  with  this;  in 
addition  the  predominant  wavelength,  which  occurs  under 
these  sea  conditions  at  an  encounter  frequency  of  -1.6 
rad/sec,  may  be  observed  in  the  FFT. 

The  linear  system  and  observed  data  do  not  agree 
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well  in  pitch  rate,  the  linear  system  predicting  an 
increasing  response  from  the  low-frequency  end  until 
resonance,  while  observed  data  shows,  as  should  be  expected, 
a low-frequency  roll-off  to  a relatively  flat  aid-range 
response.  No  plausible  explanation  may  be  advanced  for  this 
variance  at  this  time;  however,  considerable  difficulty  has 
been  observed  by  both  Booth  and  this  author  in  correctly 
interpreting  following  seas  results. 

The  observed  pitch  rate  FFT  shows,  as  have  the 
runs  in  the  ether  two  sea  conditions,  the  mysterious  double 
peak  in  pitch  rate  resonance. 


D.  SUaMABY  CF  FfiEQaSSCY  DOMAIN  ANALYSIS 


Comparisen  cf  measured  and  linear  system  theoretical 
frequency  response  shows  significant  correlation  under 
small-sea  conditions;  nulls  predicted  by  the  linear  system 
may  be  readily  observed  in  test  craft  data,  and  the 
frequency  domain  is  relatively  uncluttered  by  noise.  The 
bulk  of  noise  in  the  system  is  not  "white”  noise,  but  tends 
to  cluster  about  sharply-defined  frequencies.  The 
techniques  developed  in  this  chapter  show  the  value  of  the 
frequency  dcoain  as  an  interesting  alternative  to  the  time 
domain  as  a means  of  model  verification  and  validation. 


figure  35  - REPHES ENTATI V3  wIDE-BAND  NOISE  SPECTRUM  Of  THE 


XR-3 
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Figure  36  - OBSERVED  PITCH  FREQUENCY  SPECTRUM, 
15  KT3,  AHEAD  SEAS 
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Figure  45  - OBSERVED  PITCH  RATE  FfiEQUENCY  SPECTRUM, 
15  KTS,  FOLLOWING  SEAS 
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Figure  47  - LINEAR  SifSTEM  PITCH  RATE  FREQUENCY  SPECTRUM, 

15  KTS,  FOLLOWING  SEAS 
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VII.  SUMMARY 


A.  CONCLUSIONS  AND  OBSERVATIONS 

The  linear  two-degree-of-freedoo  model  developed  in  this 
work  has  been  shown  to  be  a reasonable  approximation  to  both 
the  Loads  and  Morions  program  and  to  actual  observed  craft 
behavior,  in  both  the  time  and  frequency  domain;  its 
relatively  uncluttered  simplicity  has  disclose!  some 
interesting  aspects  of  testcraft  behavior,  including  the 
nulls  and  nodes  that  characterize  the  craft  frequency 
response.  Furthermore,  this  simulation  is  extremely 
economical  of  CPU  time,  requiring  less  than  three  minutes  of 
CPU  time  to  generate  the  frequency  response  curves  fcr  up  to 
100  discrete  frequencies. 

The  non-linear  behavior  described  in  Chapter  V is 
somewhat  supported  by  observations  of  testcraft  data  in  rhe 
frequency  domain  in  Chapter  VI.  This  interesting  behavior 
is  apparently  akin  to  that  of  a relaxation  oscillator. 

Noise,  which  has  interfered  with  varying  degrees  of 
severity  in  almost  all  transient  analysis  done  with  the 
recorded  testcraft  data  at  NPS,  is  reduced  to  little 
significance  using  the  frequency  domain  method  applied  for 
testcraft  data  reduction  in  this  work.  The  major  source  of 
noise  in  the  XR-3  has  been  identified  as  primarily  ignition 
noise  from  the  eight  gasoline  engines  aboard  the  craft, 
transmitted  electrically  fo  the  data  acquisition  system. 
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B.  RECOaHENEATIONS  FOR  FOIURE  STUDY 


The  linear  heave-pitch  aodel  developed  in  this  work  can 
readily  be  expanded  to  include  roll,  yaw  and  speed  effects. 
Completicn  of  this  expansion  would  result  in  an  analytical 
six-degree-cf-freedom  model,  and  lead  to  a more  thorough 
understanding  of  CAB  craft  behavior.  An  example  of  the 
potential  of  this  approach  may  be  seen  in  ref  [9J » in  which 
a rudimentary  pitch  control  system,  based  upon  this  model 
and  utilizing  optimal  Ricatti  solutions,  promises  to  reduce 
sea-state  pitch  accelerations  to  less  than  ten  percent  of 
the  open  loop  response.  The  practical  problems  involved  in 
the  implementation  of  a test-bed  ride  control  system  for  the 
XB-3  testcraft  do  not  appear  to  be  severe,  aside  from  the 
physical  installation  of  a suitable  set  of  stern-plane 
control  surfaces;  an  observer,  based  on  a linear  model, 
could  be  utilized  as  an  optimal  observer  in  a Kalman 
filtering  scheme  to  generate  some  of  the  system  states  which 
may  not  otherwise  be  suitable  for  feedback.  Since  this 
system  is  only  piecewise  linear,  the  microprocessor  appears 
as  the  most  feasible  means  of  implementing  this  control 
system,  storing  the  optimal  gains  for  several  operating 
speeds  in  a table  look-up  or  interpolation  method.  Most  of 
the  system  development  could  be  implemented  on  a 
slightly-modified  version  of  the  Loads  and  Motions  program. 
This  avenue  has  great  future  potential  for  expanding  the 
sea-state  envelope  and  improving  the  habitability  of  CAB 
craft . 

The  linear  freguency  domain  response  simulation  should 
prove  useful  in  tne  study  of  the  non-linear  behavior 
described  in  Chapter  V.  It  is  felt  that  this  non-linearity 
can  be  modelled  in  the  complex  plane  utilizing  a describing 
function.  Fortunately,  the  complex  gains  of  the  linear 
system  ate  directly  available  from  the  linear  system  in  the 
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program  developed.  This  is  an  extremely  interesting  area  of 
future  srudy. 


In  the  area  of  parametric  analysis,  such  a linear  model 
may  be  utilized  with  the  technique  of  root  locus  analysis. 
An  elementary  program  was  developed  for  this  purpose  and  is 
included  at  the  end  of  the  text  with  other  programs 
developed  for  this  investigation.  However,  time  and  space 
did  not  permit  its  inclusion  into  the  main  text. 
Furthermore,  such  parametric  analysis  was  inconsistent  with 
the  main  thrust  of  this  effort.  It  is  nevertheless  hoped 
that  the  program  will  be  of  benefit  in  future  studies. 

The  XR-3  testcraft  was  not  available  for  data  collection 
throughout  the  spring  of  1977.  The  frequency  domain  data 
reduction  techniques  developed  in  this  work  could  not, 
therefore,  be  utilized  in  anything  approaching  an  in-depth 
study  of  craft  behavior.  One  area  in  which  these  frequency 
domain  techniques  may  prove  most  useful  would  be  the 
investigation  of  actual  testcraft  gains  in  heave  at  various 
frequencies,  to  ascertain  the  significance  of  wave-plenum 
chamber  interactions.  The  15  dE  discrepancy  between  the 
Loads  and  Motions  program,  which  assumes  a plenum  wave,  and 
this  linear  model,  which  does  not,  indicates  this  to  be  a 
vital  factor  in  the  accuracy  of  any  model. 

A serious  effort  at  ignition  filtering  is  certainly  in 
order.  Ignition  suppression  and  power-line  filtering  is  a 
relatively  straightforward  approach  that  can  greatly  enhance 
the  quality  of  data  collected  from  the  XR-3. 

An  additional  technique  for  data  reduction  to  the 
frequency  domain  may  become  available  in  the  very  near 
future.  The  MDS-8  microprocessor-based  data  acquisition 
system,  developed  by  Englehart  in  ref  will  be  capable 
of  sampling  sequentially  and  storing  in  mass  memory  up  to  1h 
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analog  channels,  with  a sampling  rate  up  to  7000  samples  per 
second.  Upon  completion  of  the  sampling  and  storage,  the 
MDS-8  can  perform  a variety  of  data  manipulations  upon  the 
data  file,  including  an  FFT;  the  data  file,  both  in  the  raw 
form  and  after  manipulation,  is  available  to  the  NP3  IBM 
360/67  computer  installation  via  the  CP/CMS  time-sharing 
link.  The  possibilities  of  this  system  as  a data  reduction 
system  are,  compared  with  techniques  previously  utilized  for 
data  reduction,  impressive.  However,  much  remains  to  be 
developed  in  this  area. 
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APPENDIX  A 
WIDEBAND  NOISE  DATA 

This  appendix  is  a resume  of  the  more  significant 
wide-band  ncise  spectra  observed  in  the  FFT*s  of  XH-3 
testcraft  data  collected  at  Lake  San  Antonio,  20  January 
1977. 

The  following  plots  are  EFT's  of  testcraft  data,  reduced 
as  described  in  Chapter  V,  for  a maximum  upper  frequency  of 
200  Hz,  and  a frequency  resolution  of  512  samples. 

STARBOARD  ENGINE  1/2  POWER, PORT  IDLE 

1.  Pitch 

2.  Pitch  Rate 

3.  Plenum  Pressure 

PORT  ENGINE  1/2  POWER, STARBOARD  IDLE 

4.  Pitch 

5.  Roll 

6.  Yaw  Rate 

BOTH  ENGINES  1/2  POWER 

7.  Pitch 

8.  Pitch  Rate 
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9.  Soil 


10.  Roll  Rate 

11.  Flenum  Pressure 

STaSBOARD  BULL, PORT  1/2  POWER 

12.  Pitch 
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A 2 - 5TAHB0ARD  ENGINE  1/2  POMES, PONT  IDLE 


PITCH  SATE  FREQUENCY  SPECTHUH  0-200  Hz 
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A 9 - BOTH  ENGINES  1/2  PCWEB 
ROLL  FREQUENCY  SPECTRUM  0-200  Hz 
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FBEQUESCY  COilAIil  LINEAR  SYSTEM  CSMP  III  COMPUTER  SIMULATION 


This  is  a listing  of  the  CSMP-III  comparer  simulation 
used,  with  minor  variations  in  details  such  as  freguency 
range,  tc  generate  all  freguency  response  curves  cf  the 
linear  system  used  throughout  this  work. 


The  CSMP  III  program  is  not  linked  automatically  to  the 
IMSL  Subroutine  Library  at  the  «.B.  Church  Computer  Center, 
MPS,  Monterey.  To  acheive  this  linkage  and  aocess  the 
subroutine  LECT2C,  the  following  Jon  Control  Language  (JCL) 
must  be  inserted  at  the  beginning  of  the  program  deck: 


//SYiL(8  TO  CNIT  = 33?O.V._  = ;c  =C: 
//  CC  C I IS 

//  03  ,3.A=SYS1  ..-PiL' 5 

//  CO  Ci3-5  = >t-.-<,i:jN=':Y33.4 )L.S  = 


iKC  1 »0I  5*^  = SH^  tCiiN 


Y33.C5»'?1.l:^C'‘OC 
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$t$CCMINLCLJ  SYSTEM  MCDELING  PPGGRAM  III  V1M3  TRANSLATOR  OLTPLT$$$ 
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PAEAMETHIC  EIGENVALUE  CSMP  III  SOLUTION 


This  program  provides  a root  locus  solution  to  the 
two-degr  ee-of-freedoa  linear  model.  The  reserved  cSilP 
variable  TIME  is  redesignated  as  the  parameter  of  interest 
with  the  RENAME  statement,  with  the  range  of  variation 
specified  on  the  FINTIM  card.  One  additional  variable  may  be 
stepped,  using  the  PARAMETER  statement. 

The  parametric  eigenvalue  program  also  must  access  the 
IMSL  Subroutine  Library;  see  the  Frequency  Domain  computer 
program  listing  for  the  appropriate  JCL  to  acheive  this 
linkage  at  the  W.H.  Church  Computer  Facility. 
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